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EFFECT OF FINITE SPAN ON TEEE AIBLCAD DISTRIBUTIONS FOR OSCIlLATHra WHGS 
II - METHODS OF CALCULATION AND EXAMPLES OF APPLICATION 
By Eric Relssner and John E. Stevens 


SLWiAEY 


Lift and moment distributions are calculated for oscillating wings 
of finite span on the basis of the three-dimensional theory of part I of 
this report. The results obtained are compared with the corresponding 
results of the two-dimensional theory . 

Rectangular and elliptical wings of aspect ratios 3 6 are con- 

sidered and the range of values of the reduced frequency k extends 
from 0 to about 1.5. The calculations are made for various bendirig, 
torsion, and aileron deflection fiinctions and the results are given in 
tabular and In gra,phical form. It is founa that for a given wing and 
given reduced— frequency value the effect of finite span depends appre- 
ciably on the shape of the wing deflection functions. It is also found 
that for a given wiri{;i and given deflection function the finite— span 
effect decreases as the reduced frequency increases. For wings of 
aspect ratio 3 appreciable three-dimensional effect occurs for values 
of k up to about 1.0 and for i-.'lngs of aspect ratio 6, for values of 
k up to about 0.5. 

A practical scheme of calculations is described and auxiliary ta- 
bles are given for the numerical analysis of additional examples. 
Formulas are included which allow direct incorporation of the three— 
dimensironal results in flutter determinants of the kind described in 
AAF Technical Report No. k798. Examples of flutter calculations on the 
basis of these formulas are given in an appendix to this report. Change 
of the two-dimensional into the three-dimensional air forces appears to 
be responsible for flutter— speed changes of from 10 to 20 percent for 
wings, with the possibility of larger con’ections in tail flutter. 

While in the examples analyzed so far the aerodynamic— span effect 
increases the theoretical flutter speed, the possibility of the effect 
being in the opposite direction in other examples has to be considered. 
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INTRODUCTION 


The magni't'u.cLe of* the aerodynaiiiic~span. effect in ving theory 1b one 
of the prohlems for which no complete answer has yet "been fo'und. With 
attention restricted to wings of aspect ratio that is not too small (say, 
AR ^ 3) it may he said that, for wings in uniform motion, llftl^-line 
theory represents a satisfactory solution of the prohlem. For wings in 
nonuniform motion tho prohlem is considerahly more difficult and most 
attempts to analyze the aerodynaioj.c— span effect appear to he either in- 
complete or in their applicahilitj’’ restricted to special cases of the 
prohlem. 

In part I of the present report an aerodynamic theory of the oscil- 
lating wing of finite span has been given which is considered to he as 
Inclusive as lifting-line theory is for the wing in uniform motion. The 
final results of this theorj'- were formulas for the spanwise distrihution 
of air forces and moments for the four basic types of motion of fD.exihle 
wings, namely bending, torsion, aileron, and bah deflection. In these 
formulas the effect of three-dimensionality appears as a correction term 
a to the basic function C(k) of the two-dimensional theory. The 
culation of o necessitates the solution of an integral equation which 
is similar to, hut less simple than, the integral eqmtion of lilting 
line theory for the lift (and circulation) distrihution on wings in 
uniform motion. 

With a method thus established for the systematic calculations of 
finite-span corrections for the customary two-dimensional theory it 
becomes possible to arrive at statements with regard to the quantitative 
importance of the aerodynamic— span effect. It has in the past been 
held by a number of investigators that this effect, while significant 
for wings in uniform motion, is no longer so for oscillating wings. Phys- 
ical consideration of the wake pattern indicates that for a given wing 
tho actual three-dimensional flow approaches more and more nearly the 
two-dimensional pattern as the frequency of oscillation increases. It s 
thus permissible to say that for sufficiently high frequencies the 
aerodynamic-span effect is of an insignificant magnitude. This, however, 
leaves the question as to what constitutes a sufficiently high frequency 
and what is the magnitude of the effect if the frequency is not suffi- 
ciently high. It was therefore considered desirable to establish, at 
least roughly, the range of frequencies for which there is an appireciable 
aerodynamic-span effect and to indicate the nature of the effect in this 
range. The variables which are mainly involved are o,spect ratio, 
wing deflection form, and reduced frequency k = ebb/U, where (x> 
circular frequency of oscillation, b the semichord of the wing, an 
the velocity of flight. 


The calculations in this report are in part designed to permit a 
rapid estimation of the asgnitvie of the aercry-a-imic-span effeco by 
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providing for consideration the results of an appreciable nunber of sam- 
ple cases for wings of aspect ratios 3 and 6 and rectangular or elliptic 
plan form and for a range of values of the reduced frequency k. 

It is felt that a point of further interest of this work should be 
the presentation of a scheme for incorporating the three-dimensional air 
forces in flutter calculations, for those cases in which such an incor- 
poration is deemed worth while. Although the authors agree that such a 
refinement will be unnecessary* in a good many cases for which great 
accviracy is not required, they also feel that if it is desired to obtain 
theoretical flutter speeds which are within, say, 10 percent of the 
actual flutter speeds, then in most cases it will be found necessary to 
incorporate the aerodynamic-span effect in the analysis. Appreciably 
larger corrections ^d-ll often be obtained, particularly in problems of 
tail flutter. 

This investigation, conducted at the Massachusetts Institute of 
Technology, was sponsored by, and conducted with financial assistance 
from, the National Advisory'' Committee for Aeronautics. 


SYMBOLS 


U velocity of flight 

X, y Cartesian coordinates in plane of projection of wing surface 


P 

xiiy) 

xt(y) 

b 

^o 

s 

CD 

Q 




“m 


density 

coordinate of wing leading edge 
coordinate of wing trailing edge 
semi chord 

semichord at midspan 

ratio of span to chord at midspan 

circular frequency 

amplitude of three-dimensional circulation function defined by 
equations (8) and (13) 

dimensionless spanwise coordinates; y* = ” ^Ao 

dimensionJ.eas coordinate of midchord line; = (x^^ + x^)/2bQ 


k 
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k 



X 



C(k) 


0 

a 

c 

e 

d 


f 

I, m 
h 


a 


P, 7 
L 

Mp, My 


reduced $’requency; k = oib/U 
reduced frequency at midspan 

function defined by equation (10) 

amplitude of two-dimensional circulation function given by equa- 
tion ( 7 ) 

variable of integration 
function defined by equation (16) 

function defined by Theodorsen 
correction f\inction defined by equation (6) 
location of elastic axis in units of the semichord b 
location of aileron leading edge in units of b 
location of aileron hinge line in units of b 
location of tab leading edge in units of b 
location of tab hinge line in units of b 

aileron and tab overhang in units of b; I = e — c, m = f — d 

— , . io)t 

bending deflection of wings; h = h(y)e 

— / \ iot 

angle of attack of wiiogs; a = a(y)e 

« 7T icDt 

ai3.eron and tab deflection angles; p = p e 
lift per unit of span; L = L e^^^ 
moment about elastic axis per unit of span 

aileron and tab hinge moments per unit of span 


F function defined by equation ( 9 ) 

J, Y Bessel functions in customaiy notation 
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H Hanlrel functions in customary notation 

lift coefficient defined by equation (12a) 

Cp rolllng-mement coefficient defined by equation (I2b) 

P, Q functions defined by equations (32) and (33) 

0 j(y) phase angle functions (j = h, a,, P,. 7 ) 

LIFT AND MOMENT FUNCTIONS OF THE THREE-DIMENSIONAL THEORY 


According to part I of the px’esent report the following formulas 
hold for the distribution of lift L, moment M^^, aileron hinge moment 

Mp, and tab hinge moment My due to a bending deflection h, a tor- 
sional angular deflection a, an aileron angular deflection p, and a 
tab angular deflection 7 : 


ikXC(k) + Ok]'] 


2pU^bo 2 


+ rt <fi [Ik + k^a] + j”l+ik^i — a^ j [C(k) + o^^ .]'| a 


j"^o 




[IkCp^ - + 


ik _ 1. 


El + ~ Eg j[C(k) + Op J J; p — 


r« 


I [ikC^^ - + j^Ei(d) + Eg(d) I [C(k) + Oy ] |7 (1) 
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5? - ^ f - (i* 0 ] ?/ ■ ( ^ T ’’ 


2pU^h^ 


X ( ^ ( I ~ “ ) 1 -K>a] ] O' 


+ r^rAo -fi 


[Ap^ + ikAp^ - k^Ap^ ] 


to / '2 


(^"■v 

El + -- Ep 

V 2 / 

2 


[C(k) + Op ] I p + 


4| (r" 


Ei(d) + EsCd) ][C(k) 


+ o^] ^ 7 


(2) 


Mr. 




?pU 2 k 2 bg , 2 


■"o '^'O 


J~ [lkB„ - k% ] + - E, 
1 2 “■1 2 ' 


t, 1 + Ik (I - a) I [ C(k) + 0 ^] 'l a 




h 1 - Es 


[C(k) + Op] 1= 3 


*it)° 


k% ] + ^ E 3 

^1 -' 2 rt ® 


Ei(d) + E2(d) 


X 


[C(k) + 0 ^] ] 7 


(3) 
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H 


Hankel functions in customary notation 


lift coefficient defined by equation (l2a) 


rolling— mement coefficient defined by equation (I2b) 


P, Q functions defined by equations (32) and (33) 

0j(y) phase angle functions (j = h, a, 3 . y) 

LIFT AND MOMENT FUNCTIONS OF THE THESE-DIMSNSIONAL THEOEY 

According to part I of the present report the following formulas 
hold for the distribution of lift L, moment aileron hinge moment 

Mp, and tab hinge moment My due to a bending deflection h, a tor^ 
sional angular deflection a, an aileron angular deflection 3, and a 
tab angular deflection y: 



+ ilc[C(k) + Oh] ] ^ 


j ^^o 


+ TC 





( 1 ) 
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2pU^b^ ^o 2 \2 J 




X ( |[i“c2 - >'X3 - (2 * )f ^ ( 5 “ “ ) I + %] ] “ 


- )“ f I [Ap, + ] 


-(hO 

f.1 


E,- + ~ Eo 
2 


[C(k) + Op ] I 3 + (^ 


t V 


/I . 




Ei(d)" + EsCd) [C(k) 


+ ^^7] j y 


( 2 ) 


M, 


3 b r 3ck‘ 


2pxj2b2 , 2 




x(f - 1^X3 •" 5 ®3 P I'' ( 5 - “) i c '=(>'> “a] I “ 




t J . k E 3 


•TT T. 

El H Ep 

2 


[C(k) + Op]j3 


+ ( — ^ f - [B„ + iM - k% ] + i 
V b / ^ 2 ^3 ^2 >^1 ■' 2 


2n ® 


Ei(d) + E2(d) 


[C(k) + 0^] 1 7 


(3) 
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J?l_- . L. f . 3^! D, * i E fd)lk[0(k) + o.] ) i*. * 

2pU-'bg I 2 2 3 ' ^ M b_ Vb.y 


; 


X [f ~ kX, ] + I E^Cd) |l + ik ( i - 


a ^ j [C{k) + ] j a 


r rY ^ ! [%3 + + ” 2^(4) k + ^ E2 


[C(k)_+ <?p]^ 


X 3 + f ~ V r 1 r-D .,. iki) _ k^T) ] + — j:,3 

\ bo / 1^2 ^ 7a 72 n-* 2 jt 


EJd) 


Ei(d) 


+ "EgCd) [C(k) + O.A V 

P f I 


(4) 


Equations (l) to (4) employ the notation of reference 1 and the 
terms A, B, C, and D are defined there. The values of the terms E 
are 


El = Tio ~ 


Ei(d) = Tio(d) - inT2i(d) 


Es = Txi — 22Tio 


Ea((3) = Tii(d) - 2 mTio(d) [> ( 5 ) 


'Eg = Ti 2 - 2ZT20 


EsCd) = Ti2(d) - 2mT2o(d) 


and the terms T are also defined in reference 1. 

The dimensions a, c, d, e, f, I, and m (all as defined in ref- 
erence 1 and in fig. 1) are in units of the local semi chord b (b = bQ 

at midspan) and k = cob/U stands for the local reduced frequency. 
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The terms cTj(j = h, a. p, 7) which alone represent the effect of 
finite span are given by 


"J = 


C(k) + 


iJi(k) 


J^(k) - iJi(k) 


( ^(/) 


(6) 


The two-dimensional circulation functions Q in equation (6) 

are defined by 


(^) _ lnc(lt) IV h(y) 

"h -“3-— '‘'V 

kHj (It) 


n 


(2) UiCfk) r. 

a (2) 

kHi (k) 


r / 1 

\ 

1 

1 + ik ( — 

- a ^ 

) 

L ^ 2 

J 

I 1 


j^( 2 ) ^ 4 iC(k) 


3 (2) 

«klia (k) 


El + “Es]— P(y) 

2 J bo 


Q 


(2) 


^iC(k) 


(2) 

jtkHi (k) 


Ei(d) + ^ EgCd) 


~ 7(y) 

^o 


^ ( 7 ) 




In order to take into account the possibility of phase variation in 
span (ilrection the fxmction h(y), a(y)^ 3 (y), and r(y) in equations 

-iJ^uCy) -i^w(y) 

(1) to ( 4 ) and (7) must be multiplied by terms e , e , 

-i0R(y) -i0y(y) , , 

e ^ , and e ' , respectively. In equations ( 7 ) the quantity 

kjjj = k(x^ + x^)/ 2 b^ represents the effect of sweep of the midchord line 
of tlie wing. 

The three-dimensional circulation functions Qj are solutions of 
the following integral equation: 
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fl 



— s 


In equation ( 0 ) y* Indicatee the location of points alomj the span in 
units of tho semichord at midspan 

The function F in the second integral is defined as 


V r* — iX, i 1 1 

F(x) = ' , i j - 


/ 2 2 
x\ 


d X, 


( 9 ) 


The function p is defined in terms of Bessel functions as 


u(k) = 


Jo(k) - IJi(k) 


nk J j J^(k) - Yi(k) I - i I Ji(k) + Yo(k) ] "i 


( 10 ) 




J J 


Numerical values of the functions C(k). p(k) , 


C(k) + — 


ic(k) 


J^(k) - iJ.(k) ^ 4^)(k) ' 


and F(x), which are of importance for the evaluation of the general 
formulas J are given in eahles I to V. 
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Lift and moments at points of no ving deflection ,- At points of no 
vlng deflection equations (l) to (4) cannot to used directly because one 
or more of the Oj terms hecomes infinite, while the oorrespondi.ng de- 
flection is zero. In view of equations (6) and (7) there nay to written 
for a point where a deflection fimction j(h, a, y) becoiaos ?.ero, 


C(lc) 


J^(k) - iJi(lr)‘ 




2pU^to 


4 e 


31T 


lU 


Q 






- Ch^)] 


qos) .+ 


IJi(lc) 


Jprr- iJa(k) _ 


Ml 


P,J 


E, 


3 b 


2pU^t‘ 


2 t 


C(k) + 


iC(k) 

. 


iJi(k) 




o k- 


" *^q(E) — iJi(k) _ 

Ic^) 


kn[^'(k) 


J 


> (11) 






2pU^t2 


t. 


C(k) 


LOi(k) 


Jo(k) - IJi(k) _ 


,, iC(k) 

4 — - — : e 


Q 


kfii^^Ck) 


Expressions for lift and rollln(;;-moment coefficient s.— The coeffi- 
cient of wing lift may be defined as 



sbo 



Ldy 



8st' 


w 




__L 

2pU''bo 



d 


(12a) 
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For a rectangular wing 



For an elliptical wing 

® 

The wing rolling-moment coefficient may he defined as 


S- 


J ~eb„ 


L ydy 


pU^(2ebQ)S^ 


Sv 


2pU-'ho 


(--) d (I2h) 

V sbo / ^ sIq/ 


Corresponding formulas may he written down fox* a pltching--monent coeffi- 
cient and for control-surface hinge -anient coefficients In the ^reeant. 
paper, however, use is mad© only of the definitions of Cj^ and C|^, 


SOLUTION OF TEE INTEGRAL EQUATION FOR THE CIRCULATION FUNCTION 


For the solution of equation (8) a procedxixre is chosen which is 
analogous to one of the known, methods for the solution of the liftiixg— 
line equation for uniform, motion. 

Define new variables 0 and 5 by the relations. 


y = s cos 0, T)* 


s cos f 


and write 


n j(y*) 



si n n 0 


n 


(13) 
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In termB of the new variables equation (8) becomes . 


y 


sin n 0 -b . 

Knj ^ 

n b (3 


” y^K_ 4 - cos n 9 


nj' 


s(cOB 0 — cos O) 


de 


+ iko 


' l*(k^8|cos 0 —.COB 01 ) ^ cos ii0d0 


I cos 0 - COS (?| 


Qj^^(s cos 0) 


m 


Introduce the known foracula. 


cos n0d9 _ ^ sin n 0 


cos 0 — cos e sin 0 


and write equation (l4) in the form. 


Z f sin n 0 -k ' n 


rt r sin n 0 
sin 0 


ik^s p cos 0 — cos 0 


n .} |cos 0 - cos 9 
c 


FCk^sjcos 0 - cos 9|)cos n9d0 


1 

, j 


-(2) 

= (s cos 0) 


(15) 
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Novj define a set of functions by the relation, 


S^(k^8,p) = 


sin n 0 
sin 0 


-I- 



/■ 

,/ 


cos 0 - cos 9 
|co8 0 — cos 6 


F(kQs| cos 0 


cos ©I )co8 n 0d© 


(16) 


Equation (15) may then be writton 


V 

/ 


K 


nj 


J" sin 


n 


n 


Jt D_ 

s b. 


n(k) 






Equation (17) is solved approximately by satisfying it at as many suita- 
bly chosen points as there are coefficients Write 


and 


A 


nm 


sin n 0 Jt 
n s 




(18) 


cos 


( 19 ) 


Then the system of equations to be solved is 


lU 
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Note, that if the i/ing deflection is syEmetrical atout midspan, only odd 
values of n occm-, while for anti syiranstri cal deflections only even 
values of n occur. If this is taken into account, all points 0m may 

he chosen within a semi span. If the number of points which are 

selected is N, then equation (20) is a system of N simultaneous com- 
plex linear equations to be solved for the N complex quantities K^j. 

Note further that the left-hand side of equation (20) depends only on 
wing plan form and on reduced frequency, but not on the form of the 
deflection function which determines the right-hand side. 

The functions which occur in equation (l8) have been 

calculated for a number of values of n, k s, and and the results 

are contained in table VI . The functions sin n0jj^/n are also calculated 

for the same values of n and 0 and the results are given in table 

in 

VII. Figure 2 contains a series of plots of against KqS, so that 

also for intermediate values of k s the values of S are available. 

o n 

Calculation of the functions according to equation (l6) in- 

volves the evaluation of an integral containing the function F defined 
by equation (9). With tabular values of F. avai3.able, the integration 
can be carried, out graphically or by approximating F in the X’ange of 
interest of the variable x by a fimction which permits e.xplicit inte- 
gration. Both these procedures have been used. (See appendix A.) Use 
has also been made of calculations carried out in England (reference 2) 
in connection with a theory of the aerodynamic— spein effect which differs 
from the one evaluated here . 


SUMMABY OF FBOCEDURE FOE CALCULATING THEEE-DIMENSIONAL EFFECTS 


(l) Choose spanwlse stations at which integral equation (8) is to 
be satisfied. (Present range of calculations of auxiliary 
functions restricts choice tb the odx stations, 

cos = 0, 0.2, O.k, 0.6, 0.8, and 1.0 


in units of the semispan.) 
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{2) Tabulate the values of K, kj., a, if needed, of 

dj Oj f, I and' m at these stations. Take from table IV 

( 2 ) 

the values of the function 4lC/kHi at these stations. 

, X -(2) 

Calculate, according to equations (7) the values n at 
these stations. 

(3) Calculate by means Of tables II, VI, VII the values of A of 
of equation (18). ^ 

(k) Solve the system of equation (20) for the coofflcientc 

preferably by the Crout method for the solution of systems 
of simultaneous equations with complex coefficients. 

(5) Calculate the values of Q . at the selected stations 

J 

by means of eqviation (13). itake us© of table VII for factors 
sin n0^/n. 

— ( S ) 

(6) Calculate the valiies of 0. ,/q at the selected stations. 

d J 

(7) Caiculdt© teisES Oj from equation (6) making use of tfeble III, 

(8) Substitute teimis Cj in equations (l) to (k) for lift and 

moments . 


It is this procedure that has been used to obtain the results which are 
discussed in the following section of this report. 


COMPARISON OF TWO-DIMENSIONAL AND THREE-DIMENSIONAL LIFT DISTRIBUTIONS 


In order to Judge the effect of the three-dimensional, correction 
terms in the expressions for lift and moments, a number of cases have 
been analyzed numerically and the results are shown in table VIII, and 
by means of figures 3 to 21 which contain the two-dimensional as well as 
the three-dimensional distributions. 

Calculations have been made for wings of elliptical plan form and 
for wings of rectangular plan form. In both cases the aspect ratios 
chosen are 3 and 6, and a range of values of the reduced frequency k^ 
is covered. As wing deflections the following were chosen; 

(1) Translation of rigid wii;g ^ = l.'\ 

\ °o / _ 

(2) Pitching of rigid wing about midchord (a = 1, a = O) 


i6 
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(3) Bolling of rigid ving = y/a^o^' ‘ ■ '• •• 

(4) Linear antleymmetrical torsion (a = ^ ~ 

(5) Linear eynmietrical torsicjn (ct - lyl/s^o^ 

(6) Linear syirmietrical Lending ~ 

(7) Parabolic symmetrical bending = (y/s^o^ 

(8) Deflection of full-span aileron (3-=. l) 

(9) Deflection of paj£tial-span aileron _ 

(.|yl<Xsbo : P = 0; lyj >\sbQ : p = 1) 

Fitjures 3 to 21 contain the results of these ^IculationS (Some ^ 

the reseats were first obtained in reference 3.) indicated S 

integral equation for the circulation has seen satisfied aie ^ 

the figures; the number of these points determines gl^ted^ 

neous equations to be solved. The following conclusions may be stated. 

The aeroch'namic-span effect for a given aspect ratio and given wing 
deflection decreases noticeably as the values of the reduce 
Kq increase. There is substantial agreement between the results of the 

three-dimensional theory and the results of the two-dimensional theo^ 
when the value of ICf. is sufficiently large. The rate at which this 
approach between the results of the two theories takes place ^ 

aspect ratio and %d.ng deflection function. The lower the aspect rat , 

thf slower the approach. The finite-span effect IfZn 'tlln 

k which are higher for linear spanwise variation of deflection than 

they are for no spanwise variation of deflection, and higher for 
bone variation tLi for lineai. variation, "^®^^^^^^^,^fJ!^t^a- 

would be exT<ected, the present calculations should furnish 

tive inforry.tion yoich is needed in order to k-nov when to neglect .he 
aern'-ynamic-s-nan effect and when not. to neglect it. Bodily, it may 
stated then that vdien AE = 6, there is no need to cons — • « 

and when AE = 


aerod^'namic- span effect when k^ >1 


there is no needL 
k^ the 


tc consider the effect when k^ > 2- For smaller values of 
the effect should be given greater <itr!^°’"'and’^when^k/r^^ 


certainly do occur when k 


0.5 and AE = 6 


AE - 3 Corresponding conclusions will hold for values of as^ct ratio 

“t«an 3 andlor value, .ayond these ^nl- 

numerical results obtained may^ also serve . different 

tude of the effect to be expected in cases which are somewhat different 

from those analyzed. 
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It appears from the present calculations that the aerodj"namlc-epan 
effect is appreciably less for hinge-moment distributions. It is thought 
that this fact vas not kno-'-m previously. 

In connection with present resviits for partial- span control sur- 
faces, it should be said that greater acciaracy than has been obtained is 
desirabie. The discontinuous spanwise variation of deflection should be 
taken into account either by using more terms in the series for 0 or by 
using a special function ^7 ^j_gQ ■ which accommodates the discontinuities 

in a manner analogous to that in uniform-motion wing theory. 

It is noted that for the roctan^gui.ar wings the lift distribution is 
not zero at the tip as would be expected. This may be explained as fol- 
lows: While in an exact lifting-surface theory the chordwise pressure 

distribution at the tip would vanish, this condition cannot be satisfied 
in the present approximate theory. In the approximate tlieory only an 
average condition can be prescribed and the natural choice of this average 
condition is to make the circulation vanish at the tip. For not too 
large values of k this circulation condition is effectively equivalent 
to the condition of vanishing tip-lift intensity. With increasing value 
of k those two conditions become less and less similar. The same ob- 
servation can be made regarding the moment intensities at the tip for 
rectangular wings. In an exact theory the moment functions would vanish 
at the tip while the present theory is unable to ensure this. 

It must be said that to the extent that the tip conditions are not 
completely taken care of, the present theory does not fully account for 
the effect of finite span. However, even when this difficulty occ\urs, 
corrections are obtained which are in the right direction and which 
therefore may be applied. It should be emphasized that no such diffi- 
culties occur for wings with zero tip chord because for such wings lift 
and moment intensities approach zero near the tip in the two-dimensional 
theory . 

In addition to the diagrams showing spanwise lift distributions 
there are given in figiures 15, I6, and 17 curves showing the variation 
of wing lift coefficient for rigid translation and pitching, and of wing 
rolling-moment coefficients for rigid rolling and linear anti symmetrical 
torsion. These figures show again how with the increasing k^ the re- 
sults of the two-dimensional and of the three-dimensional theory approach 
each other. 


AERODYHAMIC-^SPAN COBEECTIONS IN FLUTTEK CALCULATIONS 


The use of the finite— span corrections of this report in flutter 
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calculations inny be explained by use of the procedure described in refer- 
ence U. 

In the case of combined torsion and bending aileron flutter with 
geared tab, flutter speed and frequency are determined from a determi— 
nantal equation of the form 


The terms of this determinant Involve structural components and also 
aerodynamic components which are designated in the following calcula- 
tions by the subscript A. The aerodynamic components only are con- 
sidered here. According to and using the notation of reference h, with 
a dimensionless spanwise coordinate z, defined by 


ABC 


DBF 


0 


( 21 ) 


G n I 
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Introducing for tha L, M, T, P, and Q functions their equivalent 
values in the notation of the present repoi-t, the aerodynamic determinant 
terms become: 








V7 . \ '<( r-. '*» J 
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According to eqiiations (i) to (4) the finite— span corrections appear 
as corrections to the Theodorsen function C(k) in such a manner that 
C(k) hecoiaes C(k) + By substituting C(k) + for C(k) in the 

determinant terms, the determinant terms developed using two-dimensional 
aerodynamic considerations become three-dimensional temis. With the 
three-dimensional aerod;^Taamlc series indicated by a superscript 3 there 
may be written: 
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By rewriting the finite- span correction terms and indicating 
them hy a A prefix, the finite- -span corx-ec bions for the bending- 
tors ion-aileron flutter de ceri'iiinant of reference 4 are then as 
follows j 
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These integrals for the finite-span corrections for the terms of 
the flutter determinant cannot in general he integrated explicitly. It 
is necessary in all hut very special cases to integrate these correc- 
tions graphicta3.1y or hy numerical means. 

PftTn^i ng and torsion of rectangular wing .- This case may he consid- 
ered individually as an example in vhich explicit integration is 
possible. By taking the values of and from e(juation (o), 

there may first he written 


2i ^ iJi 
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In view of the constancy of k for rectangular wings and the fact 

that the circulation functions occur as quotients only, it is permis 
sihle to normalize the two-dimensional circulations in such a way that 


( 27 ) 
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and to indicate a corresponding change in the three-dimensional circu- 
lation functions by writing 




(28) 


Introduction of equations (2?) and ( 28 ) into equation ( 26 ) gives 
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By aseuniing linear torsion and parabolic bending, that is, 
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fjz) = z, fj^(z) = z^ 


(30) 


the remaining integrations in equations (29) icay he carried out and 
the result is 





Equations (31) will he applied in appendix B of this report for 
the flutter analysis of a rectangular uniform wing. Equations (25) 
will he applied for the analysis of a tail-flutter problem. 


THE RESULTS OF BIOT AND BOEHNLEIN 


In reference 5 theoretical expressions were given for lift ^d 
moment at midspan of a rigid elliptical wing with motion consisting 
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of translation and pitching atout the midchor'i line. The results of 
reference 5 niay he conparod wjtii corresponding roaults obtained on the 
basis of the formulas of the present report. 


The relevant expressions of reference 5 niay be written 3n the fori^ 

,2 


L(o) ^ 
2pU^b^ 


TC J - ^ + iko PAE(ko) ] ^ 


4- Tt 


, il^o N - 


2 J I 


(32) 


Mg(o) 


O ■ y o 


2pa‘^b 


+ n I i Ak, 


A„_ - kt A, 


1 2 ® o tt] j 


- . ik . . ■ 

r + -/) W^o) I 


\ 


(33) 


It may he seen that the effect of three-dimensionality in these formulas 
is responsible for the change of the function C(ho) bf the two-- 
dimensional theory into the functions and which ai'e tabu- 

lated in reference 5 Poi" a range of aspect ratio and of values of k^. 

% 

Comparison with equations (l) and (2) of the present' report shows 
the following relatione between the results of reference 5 and the 
present results: 

^Ae(^o) + '^h(^o) 


^Ae(^o) <"> C(^o) + ^^h(^o) 




(31^) 


Q;^(ko) 


<— > 


C(ko) 


+ 


Calculations show, that for the uniform deflection functions which 
are considered here, 
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ajj(ko) = ajkj,) (appro::.) 


V } 


As a formal difff^'rence between the two sets of results there re-- 
mains the fact that in the present work the correction terms for lift 
and moment are the same^ while in the work of reference 5 this is not 
so. The authors of the present paper are not able to decide which one 
of the two sets of results is nearer to the actual facts in this 
respect. It may, however, be noted that the differences between the 

P— and Q— functions are quite small. The order of imagnitude 
of the differences between the two appears to be no greater than the 
order of magnitude of the deviations of the approximate results from 
xhe results of an as yet not existing exact solution of the problem. 

A numerical comparison between the P— and Q—functions and the 
function C + a has been made for the four values of aspect ratio, 

6, 3, 2, and 1.5. Figure 22 reproduces the real and imaginary parts 
of the various functions, together with the function C of the two- 
dimensional theory. It is seen that agreement between the two kinds 
of three-dimensional corrections is quite close. A theoretically 
important difference is the fact that for increasing k^ the results 

of the present work converge toward the result of the two— dimensiona-L 
theory while, for aspect ratios 3; 2, and 1.5; the results of 
reference 5 not converge toward the result of the two-dimensional 
theorj^. As k^ approaches zero the results of the present work ap- 
proach the results of lifting— line theory. It is noted that the same 
is very nearly true for the function Q, while the function P is 
consistently somewhat smaller than would follow from lifting— line 
theory . 

Determination of terms cj .— It was found that, for the two forme 

of motion considered previously, sufficiently accurate results were 
obtained from a one— tem approximation for the circulation function 
^ . By taking 


and determining the coefficient by satisfying the integral equa— 


Q = Ki sin ^ 


(36) 



WACA TN No. 1195 


31 


and hence 


Q (o) ^ 

_ (2) 

Q io) 


1 + - |i(ko)Si ^ 




f) 


(37) 


Then, according to eouation (6) there follows for the correction teras 




- 1 


(38) 


L 1 + ^i(k^)Si ( k^s, 


Equation ( 38 ) is evaluated ty means of tables II, III^ and IV, Numer- 
ical values of the functions o so obtained are listed in table IX. It 
may be noted that equation ( 38 ), in addition to givin^^ :the correction 
term for translation and pitching^ also gives the proper correction for 
the deflection of a full-span aileron on an elliptical wing as long as 
the dimensionless aileron-hinge-line coordinate e and loading-edge 
coordinate c arc constant along the span. 


As a conclusion to this di.scusslon it may be stated that from a 
practical point of view there seems, to be little to choose oe tween the 
results of reference 5 the present results in those cases where the 
results of reference 5 apply. A practical disadvantage of the results 
of reference 5 appears to be that no provision is made to obtain corre- 
sponding' results when other than unifomm deflection functions are to oe 
taken into consideration. 

The present authors would have liked to include also a comparison 
of their results with the results of W. P. Jones in reference 2j how- 
ever, notational difficulties could not be overcome in time . 


Massachusetts Institute of Technology^ 
Cambridge , Mass . ^ January 4, 1946 . 
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APPEI3DIX A 


NOTES ON TEE E\^ALUATION OF TEE FUNCTION 
According to equation (l6) 

n 

/ cos 0 — cos d 

i 

0 

X F(kQ3|cos0— cos 0j )cos n3d0 


j cos 0 — cos 0 \ 


Sn(V> = 


sin n 
sin i 




nnd according to equation (9) 


OD 


F(x) = / 

./ 

0 



^ j dX = R(x) - il(x) 

xX ! 


The function F is given numerically in table V. It can be shown that 
it possesses the following limting behavior: 


X 



X 0 : F(x) - — In X - 0.39 Iri 2 - — 

p 


It can further be shown (reference 6) that the following representation 
holds ; 

CC 

F(x.) - i j du (a) 

• / u^ 

X 


where G(u) satisfies the differential equation, 

d^G _ 1 ^ _ Q _ i 

du^ u du u 




with the initial conditions, 
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G(0) = 1, G* (0) = -1 


(c) 


Equations (a) to (c) are listed. hers for future reference. So far, it 
has not been found possible to derive computational advantages from 
them. 

Graphical integrati on . — The integral, 

Tt 

,. 9 ,. ?^^ . , — l(k^sjcos 0 — cos 0 \ )cos n^dS 
I cos 0 - C08 0j 

0 

may readily he evaluated hy graphical or numerical (Simpson* s rule) 
me -ns as the integrand remains finite throughout. 



The integral, 
n 



cos 0 — cos G 
Icos 0 — cos 0 ! 


0 


K(kQsicos 0 ~ cos 9\ )cos n'^iO 


is subject to the difficulty that the function E is logarithmically 
infinite when cos 0 = cos 9‘, There may he written, for t 0 


7T 



0 


0 -€ 

/■ ^ 

0 


0 +e 

/" 

/ 

/ 

<p-e 




and the first and third integrals evai.uated graphically. When 0 = 0, 
write instead 



The second integral may, for sufficiently small e , be estimated 
analytically as follows: 


R(k^sj cos 0 - cos 6| ) ~ In [ 2°’ ^®k^s| cos 0 - cos e| j 
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0 = 0 + Ti d0 = cLt) 


cos 0 — cos 9 z T] sin 0 


cos n0 ~ cos n0 — ■qn sin n0 




(d) 


Hence 


0+c 

/; L 


In ( 2 sin 0| t)| ) ( cos n0 

hi 

€ 


— Tjn sin n0 ) dT| = — 2n sin n0 


r, ln( 2^ * ^ k^sr, sin 0) dri 


= — e'"^(n sin n0) j In e + ln(2° ‘ ^ '■^k^s sin 0) - 


This approximation is. satisfactory' when sin 0 kQs: < O. 05 . 

In this manner the functicns Sn have been calculated for 
kQS = 0.5 and various check calculations have been carried cut for the 

values of S obtained from reference 2. 
n 


Approximate analytical i nteg ration.— Write Sj^ in the form 


R _ = l!®!. (j/ ; i _ ^ F(ixj ) J cos n0d9 

« -/ u! j 

0 


"n 


(o) 


where 


X = kQs(cos 6 — cos 0) 


The range of integration in x is lx! <. 2kQS . In this range there 
may be written approximately: 
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la(^) 


N 

1 ^ t ^ 

— In X + / Aj^x 


n=o 


K 

/ 


n=o 


if) 


For a givon k^e an. arbitrarily close approximation for ^n may "be 
■be obtained by makinG 'the number of terms, N and M, large enough. 


For kpS < 2, (x < it), the following approximations have been 

used: 


R^(x) = -In X - 0.270 + O.T64x - 0.129x^ + O.Ollx^ 

7 

I (x) = 1.319 - O.T57X + 0.202X® - 0.02Cx^ 
a 


(e) 


The coefficients of these expressions were obtained by making R^ 

agree -vrlth R and 1, respectively, in four points of the interval 
0 < X < 4. 

By introducijjg equations (g) into equation (e) the following inte- 
grals occur: 
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cos 6 — cos -0 
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/ i cos 0 — cos 01 
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cos n0d6 = 2 
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(cos G — cos 0)cos n0d9 = 
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‘(cob 9 — cos 0) 
! cos 0 — cos 0 
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r sin(n — 2)0 2 sin n0 

\n(n-l)(n-2) n(n^ - 1) 


+ sin(n 2)0 ~) ^ ^ 
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r (cos 6 — CCS ^0 cos n0d^ = 
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8^ 


— (l + ^ cos^0), n = 1 
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— jr cos 0, n = 2 


n = 3 


0 ^ • • 


A somewhat less simple 
integral, 


calculation leads to a recursion formula for the 
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/ i 


coa 9 — cos 0 
cos 9 — cos 01 
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In I cos 6 — COB 0j I cos nddO 


j 


Tlae recursion formula is 


^ fn- - 2 cos 0 f ^ + (l + 

\ ny ^ ^ \ ny 


2 sin (n — 1)0 sin (n + 1)0 


n 


n - 1 


n + 1 




n = £, 3, 


With the initial conditions, 

f 1 = 2 sin 0 j 2 In I sin 0| + ln2 — ij — 2 cos 
fg = sin 20 I 2 in j sin 0j + ln2 — 1 j — 



— cos 


20f"0 - 

V 


2 / 


On the basis of these formulas a number of values of the functions 
Sn have been calculated and the results compared vith the results of 
the graphical solution. It vas forind that there was very good agreement 
for the imaginarj' parts of Sn, differences between the two values ob- 
tained by the two methods occurring only in the third or fourth decimal. 
Somewhat less satisfactory agreement was found for the real parts of 
Sn, with deviations of from, 1 to 3 percent. ivTaile for practical pur- 
poses an uncertainty of this magnitude with regard to the values of Sn 
is of no importance, it is believed that the results of the graphical 
analysis are the more reliable. It still appears to be desirable to 
set up an analytical scheme by means of which the values of S^^ may be 

calculated accurate to, say, three or four decimals. 

A transformation of Si(kQS, jt/2).— It may be shown that this fimc— 

tion can be written as a single integral v?-hen use is made of complete 
elliptical integrals. From eqtiation (l6) 
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5i ( I ^ 


2 ^ 


2ikoS I • — ikosX i jt 1 1 t 

' 0 


y 1 + x2 e/._ 

V -\i' 


1 :- X^ j 


dX 


where E designates the complete elliptical integral of the second Kind. 
It appears that this representation might serve as the starting point 
of further analytical work. It is prohaole that corresponding, less 
simple expressions majr he obtained for Sj^(kQS, 0) also when n f 1 

and 0 / ~ . 
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APPENDIX B 


APPLICATION OF THE THEOPY TO WING FLUTTEP PROBLEMS 
DISCUSSION OF TEE THEORY 

It vas shown in reference 4 that the equations of motion of a wing 
in Lending and torsion aileron flutter with a geared tab can bo written 
in the form of a determinaint as 


= 1 - 

p S F 

G H I 


= 0 


(Bl) 


In the present report the elements of the flutter determinant of refer- 
ence 4 are represented as the sum of a structural term, a 
two-dimensional aerodynamic term, and a three-dimensional correction 
term as follows: 


A = Ag + ^ + A A 
B = Bg + B^ + A B 
C — C 3 + C^ + A C 
D “ D 3 4* D^ + A D 
E — Eg + + A L 
F = Fg + F^ + A F 
G = Gg + + A G 
§ = 1s + 5a + ae 
I=Ig+lA + AI 






f 


J 


(B2) 


In these terras the subscript s refers to the structural terms j 
the subscript A refers to the aerodynamic terms calculated by conven- 
tional two-dimensional methods integrated along the span* and the prefix 
A ■ refers to the finite— span correction to the aerodj'namic parts of the 
elements of the determinant. 


fpll oil Pll Mil Mil Cill Mil Hll 
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The aerodynamic parts of equation (B2) are given hy equations (23) 
and (25) of this report. The structural terms are listed helov in the 
notation of reference 4. 


A / \ 1 , / “a \ I 

A 3 = (u)r , 1 - ( -- J (1 + xghj I ; 


•CD, 


'tt' 


'I .l.O ^ 


[fij(z)] ^ dz 


s = (^^^a)r 7~^ [fh(z)] L Mz) j dz 

(Sa)r 

s = (dxp)_^ ^ [fh(z)3 [fp(z)] dz 


. ^4 ^5 


3 ^ 

0 S. 


L-h(z)j Jxp(z) ] dz 


s = j '^yfr [ i^z) 3 dz 

^ 0 (&a)r 


3 = (.urci^) i 1 


(isf I j- 


1.0 I, 


■a 


(0 y 


_r 
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Z 2 Iq + ( 

H[rp2 + (e -a)xpj ' j 7- " ' ' 


J ° (Ia)r 
I3 + fe ■ a'tSp 


[ra(z)y- dz 


9 -) 


1 


n : [X [r^^ + (f - a)xg ] | 

Z 2 s 

^ Zl (S 3 ) 


z 4 Ig + (e - a)hSg 


J Z3 .[I5 + (f - a)bSg^3 


L fj 2 ) ] C ^ 3 ( 2)3 dz 

[f3Cz)]dz 


3 = (uxp)^ / 7^^ [±h(z) ] [‘3(z)]dz + n(hX5)^_ j 

r .. o / ^ ^ 7 „ ^3 + (o ■“ a)hS3 


r 

^4 S 5 


>(B3^ 


Z 3 (S 5 ), 


[fh(z)] [-3(2)1 dz 


U[r 3 ^ + (e - a)xp] ! j 


zi [Ip + (e - a)bbp]^ 


[fjz)3 [fp(z)]dz 


r„ . 


-' 1-^4 + (f - a)bSs 


+ n ! u. + (f - a)xp 3 I / —2 [fa.(2)3 [fp(z)3dz 

4^23 [I5 + (f - a)bS 53 ^ 

= (^rp-) i'l - (1 4 igp) I J ^ Cfp(z) fdz 

P r VcD^^y VcJOn/ ^ i (Ic) ^ 


2 i (^ 3 )^. 
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In addition to the general structural terme of equation (B3)^ the 
structural terms for the special case of the uniform rectangular v/^ing 
In bending and torsion flutter with a parabolic bending and a linear 
torsion mode are listed separately as follows. The coi’responding aero- 
dynamic terms are given in equations (23) and (3I). 


=s ~ 


5«pb, 


jb^ . X '' ' ^ ^ J 




Is = 




UnpDQ ^ 


> (B4) 


2s 


a 


4«pbQ " 


E 

=S 


a 




i" 

L 


/ 


OV, \ 3 

> — ) (1 iSa) 

^ 0) 


J 


where 

M mass of the wing per unit span 

static moment of wirj^;? per unit span about the elastic axis 

moment, of inertia of the wing per unit span about the elastic axis 

In order to appl;^’ finite— span correct'ons to flutter analysis a 
procedure should be outlined. A possible method is to analyze first 
the problem in question by usiri^-- the conventional two-dimensional aero- 
dynamic terms integrated along the span. If the resulting speed occurs 
in a range deemed critical, the finite— span corrections may be intro- 
duced. 

Since the anaiy^ais using the two-dimensional values of the aerody- 
namic I’arts of the terras of the flutter determinant gives the approxi- 
mate range of frequency parameters in which flutter may be expected, the 
range of frequency parameters vrhich must be investigated for an analysis 
using finite— span corrections is minimized. In general, the nianber of 
frequency parameters investigated by using ehree-dimensional aerody-namic 
'considerations may be limited to three values when this method is used. 
In the analysis conducted it has been found that the calculated flutter 


k2 
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speed recurs at a somewhat lower value of the freauency p-araaeter when 
the finite-span corrections are applied. This should oe kept in mind 
when conducting flutter analvsis which takes into account the effect of 
a finite span. 

The procedure outlined will be used In the tvro examples presented 
in this appendix. 


EXAMPLE I 


This example is presented primarily to illustrate the method of 
applying finite— span corrections to a simple vrlng in bending aiid torsion 
flutter. The wing selected for the analysi,- was the IT— 75 i-d-ng of refer- 
ence 7. In order to supply moi’e accurate data on the characteristics of 
the wing, it was again tested in the M.I.T. flutter laboratory with 
equipmrnt which was not available when the original tests were con- 
ducted. 

The N— 75 wing is a rectangular wing of aspect ratio 6 with the fol- 
lowing charac teri s t i c s : 

2sbQ = 60 In. 

2bg = 10 in. 

a = — 0 . 30 

M = 0.0CF6 siug/ft 

= O.OCO68 slug-ft/ft 

= 0.00059 slug-ft^'/ft 

a5Q^(measured \mcoupled ground 

frequenc;/ in torsion) = 8.9 cps 

(measured tincoupled groiuid 
frequency in bending) = 3-9 cps 

02^ = 0,06c 


0.070 
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■ ... Since .all experimental tests were, conducted at approximately stand- 
ard sea-level conditions, the analysis will be conducted usin^ the 
standard sea— level value for the density, 

Pq = 0,00237s slug/ft® 

- r • -• 

Nc data on the true flutter mode shapes were provided;, so the anal-^ 
ysis will he canried out with a linear torsion and a parabolic bonding 
mode aesumed. The tenms of the flutter determinant for this specr’al ' 
case have been integrated explicitly tod are provided in equations (23).? 
(31).? tod (34), The analysis will be performed and these integrated 
flutter tei7H3 used. 


Flutter .\nalj3is with '‘Vc-rimensional Values 

for the Aerodynamic Parts of the determinant 

When tho supplied data is incorpox''ated in the structural parts of 
the determinant terms from equation (34) .. these terms become^ 


y = 


M 




0^ X 2 ^ OJu X 2 _ 

: 1 - ( ~ . i — ; (i + l?h) i 

\ Cd ' V CD / 

a 


= 1.3331 J 1 - 0.1920 ^ (1 + igj^) 


is = 2 s = — ^ = 0.315i^ 


4«pb,. 


Eg = 


r / \ 2 


3jiio^o l 


= 0.8789 I 1 -(i^ V (1 + igj 

! Voj y 


The corresponding aerodj-namic parts of the determinant teimis P^p^, 
I.?. 7 and from equations (23) from tvro -dimensional considerations 
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ma;'- be calculated directly for several values" of the frequency parameter 
'oy the use of Cable 1 . • . ■ ' ' ' ' 

The two-dimensional aerodynamic terms for the rectangular wing with 
a lineal" torsion mode and a parabolic bending mode are: 










1a 



i (I. . 

k V2 




\ 

/ 


+ 


9, + 


X 

\ 




2i ,'l 
k V 2 


a 



X C(k) j 

35 Uhe aerodynamic terms oi the flutter determinant for this 
special case ma^r be calculated by simple formula substitution of the 
valixes of C(k)', a, and k, and since C(k) is tab’uiated as a function 
of k in table Ij the act'.iai calculation of the aerodynannc oerms is 
not -Derformed here. The determinant ■‘erms however, are tabulated in 
table X as a function of tlie frequency uaraiaeter k. 


The flutter determinant thus evaluated and' tabulated in table X 
will be solved by the method outlined in reference k. This me^od will 
not be illustrated, but the two solutions to each of the resultant cox^ 
plex ouadratic equations are tabulated as a f'unction of ohe frequency 
parameter k in the following tabl.e.: : - ; ^ ■ 
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ROOTS OF im FLUTTER EQUATION 


k 

First solution 

Second 

solution 

Vi 

(mph) 

Si 

(mph) 

G2 

0.4 

17.4 

-0.639 

29.0 

0.081 

.5 

13.7 

- .439 

25,1 

- .030 

.6 

11.3 

- .330 

22.3 

- .081 


The flutter speed hy this niethcd is determined graphicall;^'’ in figure 23, 
and is equal to 28.6 mph. 


Fli.rtter Analycis with Finite--Span Corrections 

Since the primary" purpose of the report is to illTistrato a means of 
incorporating finito—span correctione in flutter analysis, the means of 
calculating the cori’oction functions Cj(h, a, P, >') will be illustrat’d 
and these functions applied in this simple flutter problom. 

It has been shown in equation (31) foi* this special case that the 
finite— span corrections may be expressed in the following form; 
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Since the algebraic terms developed for this special example carjiot 
be developed easily for the genejral case, the correction temns will 
usually have to be determined bj graphical or numerical integration. For 
this example , then, the corrections will be determined by graphical inte- 
gration and by formula substitution In uhe relations listeu pre/iousay . 

It can be shown from equation (27) that the correction functxon Oj 
may be written for the rectangular wing as, 


C + 


iJi 

Jq - iJi 
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UT 


vhere f ^(y) is the deflection at the station vinder consideration and 

V is a function which is proportional to the circulation and he 
expressed as, 


w 



sin n0 
n 


The K* , terms of this Fourier series may be found from the set of 
ecuations, 





The relation, 


^ J 

p(2) ‘ f (y) 
J J 


holds only to the special case of the rectangiilar wing, but the general 
method of finding w is the same as that for finding Hj. 

In order to illustrate the method the correction terms will be cal- 
culated for = 0.4. Since the method is similar for any value of In, 
" 1 
the procedure will not be repeated for the other values of the frequency 

parameter investigated. The results, ho'vrever, will be tabulated for 
other values of k^. 

The parameters and tabulated functions necessary for the determi- 
nation of Cj for the selected value of the frequency parameter, 

kg = 0.4, are 


1)8 
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k = k = O.U 
c 

^=1 

■^o 

s = 6 

k E =2.4 
o 

M.(k) = p(k^) = 0.?644 -- Q.Q$Ghi 
C + ; = 0.5850 + 0.03091 

Jq — i«Ji 


The Sj^(0, k^s) functions listed in the following table are ob- 
tained from the plots of against kgS in fig’ore 2: 
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Wien these functions have been tabulated, the integral equation 
must be solved to fix the circulation distribution along the span. It 
is thought that the circulation can be approximated satisfactorily for 
a rectangular ving by the use of a foui>-term Fourier series approxima- 
tion. For this example the integral equation will be satisfied at 
y/sb^ = 0, 0.4, 0.8, and 1,0. 

The integral equation. 


K .A = f 
nj run 


1 

j m 


then becomes a set of four simultaneous linear equations to be solved 
for the coefficients. Since the modes assumed are symmetrical. 

only *dd values of n will occur. 


The Aj^ functions are equal to 


^Sim ” 


sin 


n 


n(]£)Sj^(0, ks) 


sb,. 


and when these terms are evaluated for n = 1, 3, 5, 7 at 

~ functions become, for y/sb^ = C, 


'•11 


= 1 + Q-Q964i) _ o.295i) 

6 


= 1 + (0.1384 - 0.05051) (0.230 - 0.2951) 

= 1.0169 - 0.05241 

A31 = -0.3333 + (0.1384 - 0.05051) (-0.493 + 0.3541) 
= -0.3831 + 0.07391 

= 0.2000 + (0.1384 - 0.0505i)(0.630 - 0.3371) 

= 0.2702 - 0.07851 
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A^j^ = -0.1429 + (0.1384 - 0.05051) (--0.715 + 0.3101) 

= -0.2262 + 0.07901 

for y/s'bQ ~ O.k, 

Ai2 = 0.9165 + (0.1384 - 0.05051) (0.280 - 0.3001), 

= 0 . 9 lf 01 - 0.05571 

A32 = -0.1099 + (0.1384 - 0.05051) (-0.204 + 0.1151) 

= -0.1323 + 0.02621 

A * -0.0936 + (0.1384 - 0.05051) (-0.326 + 0.1741) 

52 

= -0.1299 + 0.04051 • '> 

A^g = 0.1380 + (0.1384 - 0.05051) (0.772 - 0.3151) 

= 0.2289 - 0,08261 

for ?/sI)q = 0.8, 

A = 0.600 + (0.1384 - 0.05051) (0.588 - 0.080i) 

13 

= 0,6672 - 0.06841 

A = 0.U20 + (0.1384 - 0.05051X0.985 - 0.5351) 

33 

= 0.4213 - 0.12381 

A = -0.0151 + (0.1384 - 0.05051) (-0.110 + 0.1371) 

53 

= -0.0234 + 0.02451 

a^3 = -0.1398 + (0.1384 - 0.05051) (-1.330 + 0.3601) 


= -0.3057 + 0.11701 
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and for = l.O^ 

Ai4 = (0.1384 - 0.05051) (1.665 + 1.3101) 
= 0.2966 + 0.09721 

Ag^ = (0.1364 - 0.05051) (3-115 + O.T72i) 

= 0.470:. - 0.05051 

Ag^ = (0.1384 - 0.0505i)(5. 050 + 0.4851) 
= 0.7234-0.18791 

t 

A74 = (0.1384 - C. 05051) (7. 025 + 0.3441) 
= 0.9896 - 0.30721 


With the values of A detennined* and since t-he modes in ques— 

nm 

tion are Imoim, the set of simultaneous equations. 



A 


nm 



m 


may he solved. This solution is carried out in table XI by the use of 
the Grout method of reference 3. » 


VQaen the k' , te.rms have been calculated, the fxnction v vhich 

nj 

is proportional to the circulation is established and the correction 
functions = a, h) na:.’- be calculated. The values cf v at vari— 

o 


ous stations along the span are calculated in table XII and the 
= CLj h) functions are calculatod in table XIII. 


By using the correction function, the corrections to the determi- 
nant terms mavV be evaluated graphically . In table Xn the integrands 
of the follcwii^ pertinent intOvgrals are tabulated for this graphical 
integration: 
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4 , 

dz 



3 , 

cr z dz 
a 



3 , 

a.z dz 
n 





0 

a z^dz 
a 


These integrals are evaluated graphically in figure 2h, (a), (h), (c),. 
and (d) . 


For this special case it is possible to compare the values of the 
graphical integration vrith those obtained by substituting values in the 
algebraic equations of equation ( 51 ) as is lone in the folloving table: 


For = 0.4. 



The asrod;'namic correct.ions themselves are tabulated for a range 
of in table U along vit.h the calculated determinant terms cor- 

rected for finite— span effects. 

The flutter determinants corrected for finite-span effects, the 
terms of which are listed in table Z2T , can be expanded and solved. This 
was done, and the resulting complex cuairatic equation was solved by 
use of the method of r^eference 4. These results are tabulated as a 
function of k^ in the fo-llov/ing table: 
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ROOTS OF THE FLUTTER EQUATION 


^o 

First 

solution 

Second solution 

Vx 

(mph) 

Si 

Vs 

(mph) 

Ss 

0.333 

20.8 

-0.565 

33.3 

0.035 

.40 

17.2 

-.428 

29.7 

-.058 

.50 

13.5 

-.316 

25.8 

-.120 


The actual determination of the calculated flutter speed is carried out 
graphically in figure 23 and the flutter speed is equal to 3^.2 mph. 

Dis cu ssion of the results of example I .— VQien the N— 75 wing of ref- 
erence 7 with a 10— inch chord was tested for flutter, the wing was 
observed to flutter at about 3^ niph for low angles of attack. The same 
wing was again tested by the staff of the M.I.T, flutter laboratory, 
when this wing was selected for analysis in the present report, to pro- 
vide the measured uncoupled ground frequencies in bending and torsion 
and to determine the decay curves of the uncoupled bending and torsion 
modes so the damping coefficients could be calculated. In addition, 
the wing was again tested in flutter and was observed to flutter at 
3^.2 mph. 


Since the calculated value of the flutter speed with finite— span 
corrections was found to be 3^>2 mph for a value of the frequency 
parameter of about 0.3, the check between the theoretical and 

measured values Is reniarkable. It should be noted, however, that this 
close check was entirely ’unexpected since the flutter modes used in 
the analysis were not found by conducting a ground— vibration— mode 
analysis, nor were they observed in actual flutter tests. It is be- 
lieved, however, that the modes used were reasonable. 

As an ad’terthought, the M.I.T. flutter laboratory installed end 
plates in the wind tunnel as near the tips of the wing model as was 
practicable (the gap was about l/l6 in.) in order to approximate two— 
dimensioxxal flow conditions as nearly as possible. When this was done, 
the measured flutter speed was 30-1 mph. The calculated two-dimensional 
value was 28.6 mph for a value of k^ of aboub 0.4l. This too is a 

remarkable check, considering that a truly two-dimensional flow was not 
possible under the test conditions because a slight gap was left be- 
tween the wing tips and the end plates. From this analysis it seems 
that the extra work of a three-dimensional, analysis is Justified. 
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EXAMPLE II 


In this example the tail— flutter problem of Example Wo. 7 of ref~ 
erence 4 will be analj'zed with the benefit of finite— span corrections. 
The tail in question flutters in fuselage torsion and fuselage bending. 
For the sake of the analysis the surfaces themselves are assxmied to be 
rigid. 

When this example was analj'^zed in reference h, a method of equiva- 
lent aerodynamic chords was used to calculate the aerodynamic 
coefficients. This method is unsuitable for more exact analyses and is 
replaced in this problem by a two-dimensional method of calculating the 
aerodynamic coefficients which depends on the plan form of the surfaces 
to be analyzed and in which the two-dimensional aerodynamic coefficients 
are integrated along the span. In this problem the structural parts of 
the determinant terms used will be those of Example Wo, 7 (reference 4). 

The vertical find horizontal tail surfaces used in this problem are 
shown in figures 25 and 2.6. Various other physical characteristics of 
the surfaces are listed in tables X^LC and XVTI. The values of the \m— 
coupled frequencies to be used in the analysis are 

a>Q^(assumed luicoupled fuselage side-bending frequency) = 46.1 radians/sec 


oij^(assumed vincoupled fuselage torsion frequenc 5 '-) = 65.0 radians/sec 


CDp( free— rudder freqiiency) = 0 radian/ sec 


S = Gh = 3a 


= 0.038 
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From Example No. 7 of reference 4, the structural terms used are 


^ = 

/CD \2 ! 

67.27 jl- J 1.988(1 + igjj) I ft 

II 

100.1 ft 

II 

o! 

1.731 ft 

D = 
— s 

100.1 ft 

E = 

— B 

<■ / CD- 1.2 

1891 ! 1 - ( — ) (1 + iga) ! ft 

L \ CD / J 

F = 
-8 

23.86 ft 

G = 
— s 

1.737 ft 

H = 
—8 

23.916 ft 

I = 
—8 

5.992 ft 


The fuselage torsion and side— bending modes assumed in this example 
give the tail surfaces the following modes: The horizontal tail will 

flutter with a rigid rolling mode. The vertical tail will have a rigid 
rolling mode and a uniform pitching mode about the fuselage bending 
axis. The inadder will have a xmiforta deflection mode and since a geared 
tab is provided on the rudder, the tab will also have a imiform deflec- 
tion mode. It is for these modes that the aeroid^namic parts of the 
flutter determinaixt must be calculated. Beca’ase of the difference in 
the root chords, the value of jJq for the vertical tail will be l.l6 

times the of the horizontal tail. 
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Flutter Analysis with Two-Dimensional Values 
for the Aerodynamic Paxts of the Determinant 

By rearranging the relations for the two-dimensional aerod;>''namic 
terms in equation (23) to fit the special conditions of this example 
(where the H subscript indicates values pertaining to the horizontal 
tail and the V subscripts^ to the vertical tail), the terms become 


A. = (sb^) / ( 

-A o vJo \ 3.33 
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(sb ) / 

° Y J o 


3.33 ^ ^ L kf k^ Pi J 


E3 I 2Ei 


2Jt- 


'2 L 


^ V - 


C(ky) p dz 


by r r B73 iby 
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Tliese integrals may be evaluated graphically for a given value of 
the frequency parameter. By noting that the frequency parameters are 
different for the vertical and horizontal tail surfaces, the function 
C(k) can be foxind from the tables. The functions A, B, and C can 
be found in reference 1. The function E is defined in equation (5). 
There is, then, no obstacle to this graphical integration. 

The results of evaluating these integrals for a range of k^ are 

tabulated in table XVIII along with the elements of the determinant. 

The determinants themselves are solved by the method of reference k and 
are tabulated in the following table: 


BOOTS OF FLTJTTEE BETERMINAIIT FOE VAEIOUS k^ 


k XT 

oE 

^ov 

H 

1 ^ 

1 


V£ ! 

(mph) j 

j 

0.686 

0.795 

172 

0.004 

i 

268 1 

1 -0.176 1 
1 

.600 

.696 

194 

.011 

314 

-.248 1 

1 

.436 

0 

0 \ 

257 

CM 

\T\ 

0 

431 

-.303 i 

1 1 
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The calculated flutter speed is determined graphically in figure 2T. 
The true airspeed calculated by this method is 243 mph. 

Flutter Analysis with Aerodynamic Terms Corrected 
for Finite-Span Effects 

The relations for the finite-span corrections from equation (25) 
modified for this examnle become 
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A B = (sb ) / 


A C = (ah) 


o V 


+ n(sb^) 


1.0 

/ 

(■’ Yr 

c.rr ' 

Jo 

\ 3.33 / L 

kjj J 

1.0 

V 

y;'- 

3.33 i 

i + 2 i j 

1.0 / 

V -.3 r 

\ i 

3.33 / ' ■ 

1 r 2Ki 

,-.24 .. 
/ ( 

' ^ fJ\ 

1 r 2 Ei(<i) 

'h V 

. 3.33 y y 



z^dz 


a I 


”a!f j 


ICa y zdz 

J ''V,r 


+ — ; ; ( zdz 

% J 


A D = ( sb ) / 

- °Yh 


rl.D / 3 ^ 

^ •« i / 1 ^ ^ 2 i „ 

I "i ( — + a ; — a, „ zdz 

3.33 / V V 2 / k„ 4/ j 




A E = (sbo) 


\ Y jY 1 „ ; 2 2i { 1 

I — r- ) i 1 ~ + a ; + — { a ; 

) \.3.33/ i\2 ■ i- icY 


% 
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A 


p.c , 


\ 4 
j 

C / r 

■ ^1 ■■ ■ i \ •> 

L ^ ; 

J • 

J 

^ 1 — a i — ft 

V \2 / \ Jt / L 


4 




dz 


\ /i ^/ i'\f 2Ei(d) tEsCd) - 


+ n(sb ) j ( ) ■< ^ i + a ) { - ) 

V -^o ^ 3.33 / I \ 2 ./ V rt /• 




1 

X Op )■ dz 
py J 


A a = (3bo)„ ] 


r 1- 
o 


/ by I *3 i 1 

■— ■ ■ 1 — ojiY za.z 

M.33 -^ rt J 


A H = (sb ) / 


r 1 . o /• 


i V 


\4 


° Y ^0 ^ 3.33 


/ 


.. / 


It L ;.2 ir 


+ — I a ■" ds 

V ^ o • 

A/- \ / j 


ri.o / V 4 g - 22 - ■', 

A I = (sb ) j ( — - ■■ • ■; — — + : I Op /' dz 

V c k j • _) J •'■ V L Ky. Y > 


1 . •:• 


+ n(sb^) / 

V ^ o 


^Y 


4 <- 

; ^ 


3.33 


S3 rsSiCd) iSg(d) 

— I + 


2rt 


kir 


' a„ r az 
J '\T ■' 


Vdien these terms are investigated, it is observed that five types 
of Oj appear in the different relations. This means that five sets 


of simultaneous equations must be solved to determine the circulation 
distributions over these surfaces and thus the correction f\mction8 Oj 

for each value of o>/u selected. 


For the siu-faces provided in this problem, it is felt that a three- 
term Fourier series approximation for the circulation will be satisfac- 
torj'. This means that five sets of three simultaneous equations must be 
solved. 
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In this problem there are three basic t"pes of solution in each of 
which the terms are the same: the rolling of the horizontal tail 

surfaces, the rol.ling of the vertical tai.l surfaces, and the pitching 
of the vertical tail surfaces, rudder deflection, and tab deflection. 

Of the five solutions to be conducted, only the solution for the tab de- 
flection is problematical, because of the discontinuous angles of attack 
at the tab extremities which make it impossible to obtain a satisfactory 
circulation distribution with a finite number of terms in the Fourier 
series approximation. The stand is taken that some correction is better 
than none "though: so the three— term approximation for the circulation is 
used for this mode too. 

For the two rolling modes the integral equation will be satisfied 
at y/sb^ = 0.4, 0.6, and 0.3, and for the three symmetrical modes 

,v/sDq = 0, 0.4, and 0.8. For the rolling modes, n = 2, 4, 6, and for 

the sjunmotrlcal modes n = 1, 3, 5- 

In the analysis the vertical tail is handled as if It were half of 

a wing "i-mth a full-span aileron. This is assimiod to be satisfactory 

since the horizontal tail acts as an end plate at the root chord. 

In order to illustrate the procedure a value of cd/U will be se- 
lected which will give for the vertical tail a value of = C.p06 and 

for the horizontal "oail a value of k^ = 0.436. 

For the rolling of the horizontal tail the correction functions 

0 , mav be calculated as follows: 
h 

The necessary parameters and functions for the solutions are first 
listed in the following tables: 


(k^ = 0.436) 


i 

.l! 

1 

I k 


iC(k) 

iJi(k) 

n(k) 


^0 

ksi^^k) 

jQ(k)-iJi(k)j 

0.4 

1 

0.375 

0.860 

0.9126 - 0.l49oi 

0.598a + 0 . 0154 i 

0.2463 - 0 . 0673 i 

.6 

.336 

.770 

.9451 - .17471 

.6205 - .OOSoi 

.2834 - .10381 

• 8 | 

1 

.286 

i 

.655 

. 9953 - . 2054 i 

.6523 - . 039^1 

.3007 - .10891 
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Ir s = 1.101, j/bI = O.i*. 
o ' o 


n 

sin n0 
n 


2 
1 4 
6 

0.3667 
-.2493 
.1132 ' 

i ' 

0.494 - 0.273i 
-.850 + .2931 
i .575 - .150i 

j/s^o = 0.6 

1 

2 

h 

,..ij 

0.4300 1 0.780 - 0.390i 

-.1344 i -.547 + .I60i 

1 -.1098 ; -.710 + .176i 

1 1 

1 

! = 0.8 
1 

j 

n 

2 

k 

i 6 

1 i 

0 . 4800 

.1344 

-.1098 

L.. 

0.768 - 0.4331 i 
.525 - .1701 I 

-.710+ .20311 

1 


V/hen this is done, the coefficients of the integral equation, 


\ 


\ 





5 m 


ma^' be evaluated. 

3y taking the necessary values from the tables and noting that 


A 


nm 


sin n0 


+ 


n 


“ h(k)Sj^(0, k^s) 

eb^ 


the Ajy^ coefficients become. 
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for y/sb^ =0.4, 

Api = 0.3667 + (0.2453 - 0.08731) (0.494 - 0.2781) 

21 2.525 

= 0.3667 + (0.2635 - 0.09341) (0.494 - 0.2781) 

= 0.4709 - 0.11941 

A41 = -0.2493 + (0.2635 - 0.09341) (-0.850 + 0.2931) 

= -0.4459 + 0. 15661 

■ A_ = 0.1132 + (0.2635 - 0.09341) (0.575 - 0.1501) 

= 0.2507 ~ 0.09321 

for y/eb^ = 0.6, 

Ap2 = 0.4800 + (0.2834 - 0.10381) (0.790 - 0.3901) 

2.525 

= 0.4300 + (0.2715 - 0.09941) (0.790 - 0.3901) 

= 0.6559 - 0.18441 

A^2 = -^il344+ (0.2715 - 0.09941) (-0.547 + 0.l60i) 

= -0.2670 + 0.09781 

Ag2 = -0.1098 + (0,2715 - 0.09941) (-0.710 + 0.1781) 

= -0.2849 + 0.11111 
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for j/sh^ =0,8, 

A = O.U'SOO + (0.300T - 0.10891) (0.788 - 0.4331) 

2.525 

= 0,4300 + (0.2450 - 0.08871) (0.733 - 0.4331) 

= 0,6347 - 0.1761 

A.o = 0.1344 + (0.2450 - 0.08871) (0.525 “ 0.1701) 

Tr 3 

= 0.2479 - 0,03821 

A^^3 = -0.1093 + (0.2450 - 0.08871) (-0.710 + 0.2031) 

= -0.2657 + 0.11271 

V?hen the coefficients have been calculated, the tvo- 

dimensional circulation functions must he evaluated. For the 

rolling mode, 


ikm 

o ik fj^(j') 


where in this case. 


-(a) , 

=4 
n 


10 (k) 

i (2)^ ^ 

; k El k 


l*v(j') 


y 

^^o 


In this example the sveepbaclc function e appears as there is 
some sweep in the surfaces. The trailing and leading-edge coordinates 
of the surfaces can he found in figures 25 and 26. By using the pre- 
viously tabulated parameters and functions and introducj.ng the sweep 

functions, the values of may he determined as follows: 
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for y/sb^ = O.U, 


^ = 


ikm 

e 



+ x-t _ -0.733 -t- 0.968 

2bo ' 2 


0.128 


k z = 0.0558 

o m 


(0.9985 + 0.05581) 

4(0.9126 - 0.14961) (0.9985 + 0.0558i)( 0.37501)(0.4) 
0.0591 + o. 55181 


for y/sb^ =0.6, 

, J±.j00 r „o,gJ2 , 0.1s; 

HI 2 

= O.O035 

(0.9965 + 0.8341) 

4(0.9451 - 0.17471) (0.9965 + 0.08341) (;i 0.336)(0.-S) 
0.0768 + 0.77091 

for y/sb^^ = 0.8, 

, . . 0.190 

in 2 

V = 0.0961 

= (0;9954 + 0.09601) 

= 4(0.9953 - 0.20541) ( 0 . 995 ^ + 0.0960i)(0.286i)(0.8) 
= 0.0996 + 0,92451 


Lkm 
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The eefc of equations, 


I 


K A 
nj nm 



.(2) 

hm 


may now be solved. This is done In table XIX by the Crout method. 
When the values of have been determined, the values of 



r — ) 



sin n0 
n 


can be found and this is done in table XX. 

The values of the correction functions can now be calculated with 
it noted that 



C(k) + 


iJi(k) 

J^(k) - iJi(k) 



This is done as follows. 


for y/^ib^ = 0.4, 


^hH 


(0.5988 + 0.0154i) I 


-0.0594 + 0.43421 

0.0591 + 0.55181 


= -0.1427 + O.llOi 


for y/sb^ = 0.6 


*^hH (0.6205 - 0.0o88i) 


-0.0895 + 0.5804i 
0.0768 + 0.77091 
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for y/ sIq =0.8 


= (0.6523 - 0.039«) 


-0.1057 + 0.6686i 
0.0996 + 0.92451 


An exactly siiailar procedure must be followed to find the correc- 
tions, cTp.^, and for a civen value of It shoxad 

be remembered though that the value of k^ is different for the ver- 
tical tail. 


When 1 the correction functions have been evaluated, the values 
of the finite— span corrections can be found by a graphical or numerical 
integration of the integrals listed in the beginning of this section of 
Example II. This integration was carried out, and the values of the 
finite— span corrections found are listed in table XXI. The flutter de- 
terminant terms calculated with use of two-dimensional aerodynamic 
considerations and those corrected for finite— span effects are also 
listed in this table. 


The determiriants can now be solved. For this problem the solutions 
are found by the method of reference 4 and the results are tabulated in 
the follo-;vLng table; 


^oH 


First 

solution 

Second 

solution 

Vi 

(mph) 

Sx 

Vs 

(mph) 

62 

0.436 

0.506 

267 

0.030 

423 

-0.212 

.600 

.656 

196 

.012 

308 

-.151 

.686 

.19^ 

173 

.006 

233 

-.071 


The actual flutter speed is determined graphically in figure 27. The 
calculated value of the true airspeed at vhich flutter occurs when the 
finite— span corrections are used is 299 mph. 

Discussion of Example II .— Example No. 7 of reference 4 indicates 
that the observed flutter speed for the tail analyzed in this problem 
fluttered experimentally at about 262 mph true airspeed. The fact that 
the theoretical analysis, however, does not check the observed flutter 
speed is not considered significant because the observed flutter speed 
was reported by a pilot some time ago and may not be accurate and the 
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analysis is approximate in that no attempt was made to take into account 
the effective inertia due to wing motion or the elastic restraint and 
damping in the rudder control system. 

It is considered significant that the three theoretical analyses 
give such widely varying values foi* the calculated flutter speed. These 
values of the true airspeed ex3 220 mph for flutter speed calculated by 
the method of equivalent chords of reference h; 2^3 mph for a value of 
kQ for .the flutter speed, of about 0.55, when the two-dimensional aero- 
dynamic terms are integrated along the span; and 299 mph for a value ox 
k^ of about 0.h5, when the aerodynamic terms are corrected for the ef- 
fects of a finite spein. 
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TABLE I TABLE II 

VALUES OP FUNCTION C = F + iG VALUES OP, FUNCTION y^(k) 



k. 


0.00 

. 0.5000 - 0.00001 

0.02' 

O'.kSio - 0.0I1221 

0.0k 

o.k6o7 - 0.06651 

0.06 - 

0-.1410 - 0.08291 

0.08 

0.1+226 - 0.09li2i 

0.10 

• o-.k05i - o.ioi8i 

0.125 

0.5857 - 0.10821 

0.15 • 

' . 0.5690 - 0.11201 

0.175 

0.5522 - 0.11501 

0.20 - 

0.5595 - 0.11591 

0.225 

0.5268 - 0.11521 

0.25- 

o.3i^k - 0.11161 

0.275 

0.50k9 - 0.10991 

0.50 

•0.2955 - 9.10761 

0.55 

0.2787 - 0.10251 

o.ko 

0.26M4. - 0.096ki 

0.14.5 

0.2519 - p. 09051 

0:50 > 

•• 0-.2k08 - 0.08k2i 

•0.55 ^ • 

0.2289 - 0.07751 

0.60 

'■ 0.2220 - 0.07221 

0.65 

0.2159 - 0.06651 

0.70 ■ 

0-.2062 - 0.O6lOi 

• 

0.1991 - 0.05571 

0.60 - 

0.1921+ - 0.05071 

0.90 •' 

‘ 0.1801 - 0.04151 

1.00 

0.1688 - 0.05291 

1.25 

* 0.1I1.56 - 0.01591 

1.50 

0.1218 - 0.00421 

1.-75 • 

0.1027 + 0.00501 

2.00 

0.086I+ + 0.00661 

2.10 

0.0807 + 0.00701 

2.15 

0.0780 + 0.00721 

2.20 

0.0751+ + 0.00721 

2.25 

0.0750 + 0.00721 

2.50 

0.0706 + 0.00701 

2.55 

o.u68U + 0.00691 

2.k0 

0.0665 + 0.0066i 

2 .I4.8 

0.0652 + o.coc:: 

2.5I+ 

o.ooio + 0.00571 

2.60 

0.0591 + 0,00521 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS. 


NACA TN No. 1195 


TABLE III 


VALUES OP FUNCTION 

* J„(kt-ii7Tk7 


k 

C + 

* Jo -IJ^ 

0.00 
0.01 
0.025 
o.ol+o 
0.05 
0.06 
o.o8 
0.10 
0.12 
o.i6 
0.20 
0.2k 
0.50 
0,3k 
o.ko 
o.hk 
0.50 
0.56 
0.60 
0.65 
0.70 
0.75 
0.80 
0.85 
0.90 
0.95 
1.00 
1.20 
1.50 • 
2.00 • 
5.00 • 

1.0000 - O.OOOOi 
O.982I4. - 0.0M2i 

0.9545 - 0.07471 

0.9265 - 0.09601 
0.9084 - 0.10551 
0.8911 - 0.11261 
0.8588 - 0.12041 
0.8295 - 0.12241 
0.8027 - 0.12021 
0.7564 - 0.10791 
0.7176 - 0.08911 
0.6845 - 0.06711 
0.6425 - 0.05101 
0.6180 - 0.00651 
0.5850 + 0.05101 
0.5646 + 0.05541 
0.5554 0.09151 

0.5075 + 0.12591 
0.4889 + 0.14851 
0.4658 + 0.17551 
0.4424 +.0.20121 
0.4187 + 0.22621 
0.5945 + 0.24981 
0.5706 + 0.27251 

0.5445 + 0.29551 

0.5179 + 0.51551 
0.2909 + 0.55191 
0.1747 + 0.59091 
•-0.0220 + 0.42451 

•-O.555O + 0.27971 
--O.1255 - 0.52291 
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TABLE IV 

VALUES OF FUNCTION 


kH; 



iC , 

k 

k Hjf ^ 

0.00 

1.5708 - O.OOOOi 

0.01 

1.5420 - 0.07601 

0.04 

1.4515 - 0.18001 

0.05 

1.4218 - 0.20121 

0.06 

1.5952 - 0. 21881 

0.08 

1.5595 - 0.24291 

0.10 

1.2901 - 0.25671 

0.12 

1.2451 - 0.26441 

0.16 

1.1671 - 0.26511 

0.20 

1.1020 - 0.25071 

0.24 

1.047^ - 0.25191 

0.50 

0.9795 - 0.19751 

0.54 

0.9415 - 0.17221 

0.40 

0.8921 - 0.15541 

0.44 

0.8652 - 0.10751 

0.50 

0.8241 - 0.06911 

0.56 

0.7889 - 0.05191 

0.60 

0.7671 - 0.00771 

0.65 

0.7411 + 0.02151 

0.66 

0.7562 + 0.02751 

0.70 

0.7164 + 0.04991 

0.75 

0.6921 + 0.07721 

0.80 

0.6691 + 0.10541 

0.85 

0.6471 + 0.12891 

0.90 

0.6258 + 0.15281 

0.95 

0.6015 + 0.17581 

1.00 

0.5791 + 0.19781 

1.20 

0.4901 + 0.27551 

1.50 

0.5544 + 0.56o6j 

2.00 

0.1281 + 0.42091 

5.00 

- 0.2279 + 0.27991 
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VALUES OF FUNCTION F(x) 


TABLE V 



1 y/x 


A 


X A 




X 

0.00 




F 


.05 

.10 

.20 


.^0 

• I<.0 


•iS 

.90 

1.00 

1.10 

1.20 

1.50 

1.14.0 
1.50 
1.60 

1.70 

1.80 

1.90 

2.00 


2.750 - 
2.109 - 
I.I1.9O - 

1.155 - 

0.955 - 

Ui : 

0.5614. - 
o.I|.89 - 
0.U27 - 
0.576 - 
0.555 - 
0.297 - 
0.265 - 
0.258 - 
0.21I4. - 
0.19^ - 
0.176 - 
0.160 - 

0.li4-6 - 
O.I5I4. - 


1.5711 

1.4681 

1.5751 

1.21+81 

I.II16I 

1.0601 

0.9871 

0.9221 

0.8651 

0.8151 

0.7681 

0.7261 

0.6891 

0.65l|.i 

0.621+1 

0.5951 

0.5671 

0.5I+01 

0.5171 

0.4961 

0.1+751 

0.1+581 


X 

F 

2.2 

0.115 - O.i +251 

2.4 

0.097 - 0.5951 

2.6 

0.085 - 0.5691 

2.8 

0.072 - O. 5 I+ 5 I 

5.0 

0.065 - 0.5241 

5-2 

0.055 - 0.5051 

3 -k 

O.OI +8 - 0.2891 

5.6 

O.OI +5 - O. 27 I+I 

5.8 

0.059 - 0.2601 

4.0 

0.055 - 0.2481 

1+.2 

0.051 - 0.2571 

i^.J+ 

0.028 - 0.2261 

1+.6 

0.026 - 0.2171 

1+.8 

O.O 2 I+ - 0.2081 

5.0 

0.022 - 0.2001 

5.2 

0.020 - 0.1921 

5.4 

0.018 - 0.1851 

5.6 

0.017 - 0.1791 

5.8 

0.016 - 0.1721 

6.0 

0.015 - 0.1671 


l/ 2 x® l/x 


TABLE VII 
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VALUES OF FUNCTIONS AGAINST cos ^ 

n 8 


^ =COS 

0 

0.2 

0.1+ 

0.5 

0.6 

0.8 

0.9 

0 

1.00 

sln^ 

1 

1 

0.9798 

0.9165 

0.8660 

0.800 

0.600 

0.4555 

0 

sin 2 /( 

■ 2 

0 

. 0.1959 

0.5667 

0.1+550 

0.1+800 

0.4800 

0.5921 

0 

sin 5)^ 
5 

-0.5555 

-0,2714]+ 

-0.1099 

0 

0.1175 

0.5120 

0.5255 

0 

sin 4 ^ 

0 

-0.1805 

-0.2495 

-0.2165 

-0.151414- 

0.1541+ 

0.2455 

0 

sin 

5 

0.2000 

0.1070 

-0.0956 

-0.1752 

-0.199U 

-0.0151 

0.1550 

0 

sin 6^ 

" T- 

0 

0.1558 

0.1152 

0 

-0.1098 

-0.1098 

0.0707 

0 

sin 7^ 

7 

-o.ii+29 

-0.0I+72 

0.1580 

0.1258 

0.0295 

-0.1598 

-0.0022 

0 

sin 8^ 
8 

0 

-0.1249 

0.0188 

0.1085 

0.115I1. 

-0.1154 

-0.0558 

0 

sin 9^ 

9 

0.1111 

-0 . 0246 

-0.0940 

0 

0.0980 

-0.0525 

-0.0882 

0 
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TABLE VI 


VALUES OF FUNCTION 

’ 7 T 


S„(/, k s) = 


sin 


n 


sin 



2 ^ + 1 / COS^- C 03 8 „ „ g 

op jr y Icos/o- cos 0| 


FOR VARIOUS VALUES OF n, koS^ and yVs = cos )Z^ 


k s = 0 
o 


X. cosp 

0 

0.2 

0 .I4. 

0.6 

0.8 

1.0 

1 

1 

1 

1 

1 

1 

1 

2 

0 


0.8001 

1.2000 

1.6000 

2 

5 

-1 

-o.8I|.oi 

-0.5595 

0 . 14-599 

1.5601 

5 


0 


-1.0879 

-0.6722 

0.8965 

1^ 

5 

1 

0.5U61 

-0.5105 

-1.2I4.6I4. 

-0.1259 

5 

6 

0 


0.6792 

-0.825I4. 

- 1-5725 

6 

7 

-1 

L0.5575 

1.0514.2 

0.2578 

-1.6507 

7 

8 

0 


0.1644 

1.1555 

-1.5113 

8 

9 

1 

-0.2257 

-0.9252 

1.1020 

-0.7876 

9 
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TABLE VI (continued) 


NACA TN No. 1195 


koS = 1 



0 

0.2 

0.4 

0.6 

0.8 

1.0 

1. 

0 ,468 
-0.5561 

0.479 

-0.5521 

0.516 

-0.5401 

0.589 

-0.5091 

0.729 

-0.2251 

1.250 

0.7271 

2. 

0.000 

O.OOOi 


0.516 

-0.2741 

0.825 

-0.5821 

1.2204 

-0.4191 

2.068 

0.4701 

5 

-0.725 

0.5011 

-o.bi5 

0.2491 

-0.274 

0.0941 

0.551 

-o.i4oi 

1.280 

-0.5841 

5.021 

0.5511 

k 

0.000 

O.OOOi 


-0.869 

0.2821 

-0.556 

0.1551 

0.7610 

-0.1621 

4.009 

0.2491 

5 

0.820 

-0.2441 

0.451 

-0.1501 

-0.424 

0.1251 

-1.065 

0.2701 

-0.116 

0.0651 

5.005 

0.2011 

6 

0.000 

0.0001 


0.584 

-0.1411 

-0.720 

0.1671 

-1.269 

0.1891 

6.005 

0.1661 

7 

-0.868 

0.205 

-0.516 

0.0521 

0.925 

- 0 . 206 i 

0.251 

- 0 .o 421 

-1.498 

0.2261 

7.002 

0.1451 

8 

0.000 

O.OOOi 


0.146 

-0.0521. 

1.026 

-0.1871 

- 1.402 

0.1751 

8.001 

0.1251 

9 

0.895 

-0.1771 

-0.201 

0.0661 

-0.854 

0.1541 

1.008 

- 0.1711 

-0.756 

0.0671 

9.001 

0.1111 


koS = 2 



0 

0.2 

0.4 

0.6 

0.8 

1.0 

1 

0.271 

-0.5181 

0.284 

-0.5201 

0.527 

-0.5211 

0.421 

-0.5181 

0.622 

-0.2691 

1.551 

1.1751 

2 

0.000 

0.0001 


0.551 

-0.2961 

0.595 

-0.4291 

0.97? 

-0,5041 

2.210 

0.8771 

5 

-0.545 

0.5551 

-0 .468 
0.2941 

-0.221 

0.1151 

0.245 

-0.1751 

1.062 

-0.5051 

5.085 

0.6471 

4 

0.000 

0.000 


-0.704 

0.5611 

-0.472 

0.2001 

0.651 

-0,2051 

4.038 

0.4961 

5 

0.679 

-0.5211 

0.576 

-0.1721 

-0.552 
0 .1661 

-0.912 

0.5701 

-0.112 

0.1151 

5.014 

0 .4051 

6 

0.000 

0.0001 


0.505 

-0.1951 

-0.650 

0.2581 

- 1.179 

0.2821 

6.012 

0.5321 

7 

0.2881 

-0.299 

0.0441 

0.814 

- 0 . 294 i 

0.208 

-0.0581 

-1.377 

0.5511 

7.006 

0.2851 

8 

0.000 

0.0001 


0.150 

-0.0471 

0.951 

-0.2751 

-1.301 

0.2521 

8.004 

0.2501 

9 

i 

0.805 

-0.2591 

-0.178 

0.1081 

-0.755 

0.2271 

0.925 

-0.2581 

-0.687 

0.0861 

9.003 
0.2221 1 
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ICqS = 



0 

0.2 

O.k 

0.6 

0.8 

1.0 

1 

0.152 

-0.2514.1 

O.II4.O 

-0.2561 

0.16k 

-0.2k8i 

0.256 

-0.2821 

o.kBk 

-0.5151 

2.055 

1.7621 

2 

0.000 

O.OOOi 


0.195 

-0.2514.1 

0.590 

-0.5901 

0.689 

-0.5251 

2 . 5 k 0 

i.k 99 i 

5 

-0.514.1 

0.5371 

-0.500 

O.28UI 

-0.160 

0.1201 

0.159 

-0.1701 

0.766 

-0.5721 

3.27k 

1.20ki 

k 

0.000 

O.OOOi 


-O.U88 

o. 58 ki 

-0.359 

0.2521 

O.U68 

-o.i8ki 

k.ik 5 

0.9651 

5 

0.14.86 

-0.5521 

0.275 

- 0,1911 

-0.250 

o.iBki 

-0 . 688 
o.k 55 i 

-0.11k 

0.2221 

5.065 
0 . 8oki 

6 

0.000 

0.0001 


0.589 

-0.2521 

-0. 14-87 
0.2961 

0.5881 

6.ok6 

0.6581 

7 

-0.58I4. 

0.5511 

-0.273 

0.0514.1 

0.6k0 

-6.5671 

0.176 

-0.0711 

-1.172 

o.k 3 ii 

7.025 

0.5651 

8 

0.000 

0.0001 


o.ioi 

-0.0651 

0.772 

-0.5661 

-1.120 

o.5iki 

8.01k 

o.k 97 i 

9 

0.651 

-0.5371 

-0.114.0 

0.1731 

-0.625 

0.5021 

0 . 77 k 

-0.5551 

0.0881 

9.008 

o.kk-2i 


k s = 6.0 
o 



0 

0.2 

o.k 

0 . 6 

0.8 

1.0 

1 

0.082 

-0.1891 

0.086 

-0.1871 

0.109 

-0.1781 

0.179 

-0.2k01 

0.587 

-o. 55 ki 

2 .k 7 k 

2.2571 

2 

-0.000 

0.0001 


0.126 

-0.2171 

0.2k2 

-0.5351 

0.511 

-o.k 97 i 

2.862 

1.9891 

5 

-0.2k0 

0.5001 

-0.215 

0.2561 

-0.125 

0.1181 

0.089 

-o.ik6i 

0.579 

-0.5691 

5.502 

1.6721 

k 

0.000 

0.000 


-0.56k 

0.5601 

-0.28k 

0.2601 

0.541 

-0.0801 

1.5861 

5 

0.571 

-0.5521 

0.210 

-0.1851 

-0.188 

0.1691 

-0.556 

0 .k 55 i 

-0.119 

0.5591 

5.ik8 

1.1901 

6 

0.000 

0.0001 


0.508 

-0.2561 

-0 .'585 
0.5061 

- 0.928 

o.l4k8i 

6.099 

0.9751 

7 

-o.k 65 

0.5591 

- 0 . 25 k 

0.0561 

0.511 

-0.58k! 

0.15k 

-0.0751 

-1.005 

o.k62i 

7 .ok 6 

0.8591 

8 

0.000 

0.000 


0.077 

-0.0681 

o. 6 k 5 

-o.koki 

-0.959 

0.5161 

8.02k 

o. 7 koi 

9 

0.551 

-0.5671 

-0.109 

0.2251 

-0.517 

0.5351 

o. 6 k 9 

-O.kOOi 

-0.512 

0.0571 

9.022 

0.6591 
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tabij: VIII 


BACA TN No. 1195 


UNIFOEM TRABSLATION OF RECTANCULAB WING OF AE = 3 


k 

7 

c 


l(2) 

L 

l)o® 

2pU^oli /bo 

2pUaboh/bo 

0.167 

0.167 

0.167 

0.167 

0 

O.k 

0.8 

1.0 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

-0.105+0.1291 

-0.157+0.1251 

-0.285+0.1151 

-0.750+0.1051 

0.055+0.5971 

0.055+0.5971 

0.055+0.5971 

0.055+0.5971 

-0.015+0.5421 
-0.011+0.5251 
-0.006+0.2471 
-0 .0002+0 . oo4i 

0.555 

0.555 

0.555 

0.555 

0 

1.0 

0.650-0.1751 
0.650-0.1751 
0. 650-0. 17S1 
0.650-0.1751 

-0.045+0.0781 

-0.069+0.0851 

-0.196+0.0851 

-0.622+0.0101 

0.008+0.6811 

0.008+0.6811 

0.008+0.6811 

0.008+0.6811 

-0.075+0.6561 

-0.079+0.6081 

-0.079+0.4751 

-0.005+0.0291 

0.667 

0.667 

0.667 

0.667 

0 

o.k 

0.8 

1.0 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

-0.018+0.0501 

-0.056+0.0551 

-0.127+0.0091 

-0.459-0.1851 

-0.426+1.1911 

-0.426+1.1911 

-0.426+1.1911 

-0.426+1.1911 

-0.489+1.1551 

-0.495+1.1161 

-0.444+0.9271 

-0.042+0.2501 

UNIFOEM TRANSLATION OF RECTANGULAR WING OF AR = 6 

0.167 

0.167 

0.167 

0.167 

0 

O.k 

0.8 

1.0 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

-0.020+0.0591 
-0.040+0.0701 
-0.169+0. 1041 
-0.750+0.1051 

0.055+0.5971 

0.055+0.5971 

0.055+0.5971 

0.055+0.5971 

0.024+0.5841 

0.018+0.5751 

0.005+0.5141 

-0.0002+0.0041 

• • • • 

1 0 0 0 0 

0 

Si 

1,0 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

0.001+0.0271 

-0.015+0.0591 

-0.100+0.0651 

-0.622+0.0101 

0.008+0.6811 

0.008+0.6811 

0.008+0.6811 

0.008+0.6811 

-0.020+0.6811 

-0.052+0.6671 

-0.060+0.5761 

-0.002+0.0291 

0.667 

0.667 

0.667 

0.667 

0 

o.k 

0.8 

1.0 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.000+0.0071 
-0. 006+0. 0l4l 
-0.058+0.0171 
-0.459-0.1851 

-0.426+1.1911 

-0.426+1.1911 

-0.426+1.1911 

-0.426+1.1911 

-0.441+1.1951 

-0.455+1.1801 

-0.461+1.0711 

-0.042+0.2501 

1 UNIFORM PITCHING OF RECTANGULAR WING OF AE = 3 (a= 0 ) 

k 

y 

c 



L 


2 pU^Qa 

2 pU^Qa 

0 

0 

0 

0 

0 

1.0 

1 

1 

1 

1 

-0.515 

-0.548 

-0.505 

-1 

5.142 

5.142 

5.142 

5.142 

2.158 

2.049 

1.554 

0 

0.167 

0.167 

0.167 

0.167 

0 

o.k 

0.8 

1.0 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

-0.105+0.1291 

-0.157+0.1251 

-0.285+0.1151 

-0.750+0.1051 

2.427-0.1501 

2.427-0.1501 

2.427-0.1501 

2.427-0.1501 

2.065+0.2481 

1.965+0.2271 

1.501+0.1571 

0.044+0.0051 

0.555 

0.555 

0.555 

0.555 

0 

o.k 

0.8 

1.0 

0.650-0.1751 

0,650-0.1751 

0.650-0.1751 

0.650-0.1751 

-0.045+0.0781 

-0.069+0.0651 

-0.196+0.0851 

-0.622+0.0101 

2..,i^5+o.5i5i 
2. 155+0.. 5151 
2.155+0.5151 

2.155+0.5151 

1.957+0.5561 

1.875+0.5411 

1.475+0.4741 

0.175+0.0221 

0.667 

0.667 

0.667 

0.667 

0 

o.k 

0.8 

1.0 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

-0.018+0.0501 

-0.056+0.0551 

-0.127+0.0091 

-0.459-0.1851 

1.924+1.2541 

1.924+1.2541 

1.924+1.2541 

1.924+1.2541 

1.857+1.5101 

1.785+1.5011 

1.517+1.1291 

0.675+0.1781 
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UNIFORM PITCHING OF RECTANGULAR WING OF AR = 6 (a = 0 ) 


0 

0 

0 

0 

0 

O.k 

0.8 

1.0 

1 

1 

1 

1 

- 0 .i 7 i ^5 

-0.2057 
-0. 514-67 
-1 

rT 42 

3.142 

5.142 
5.142 

2T5T7 

2.502 

2.052 

0 

0.167 

0.167 

0 . 1 b 7 

0.167 

0 

o.U 

0.8 

1.0 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

-0.025+0.0591 

-0.014,1+0.0701 

-0.158+0.0991 

-0.750+0.1051 

2.427- 0.1301 

2.427- 0.1501 

2.427- Q.1501 

2.427- 0.1301 

2.554+0.0451 

2.279+0.0781 

1.905+0.1591 

0.045+0.0051 

0.555 

0.5^3 

0.555 

0.355 

0 

o.k 

0.8 

1.0 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

0.001+0.0271 

-0.015+0.0391 

-0.100+0.0651 

-0.622+0.1041 

2 . 155 +Q -5151 

2.155+0.5151 

2.155+0.5151 

2,155+0.5151 

2.121+0.4011 

2.075+0.4501 

1.765+0.4681 

0.175+0.0221 

0.667 

0.667 

0.667 

0.667 

0 

0 .14. 
0.8 
1.0 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.000+0.0071 

-0.006+0.0141 

-0.058+0.0171 

-0.459-0.1851 

1.924+1.2541 

1.924+1.2541 

1.924+1.2541 

1.924+1.2541 

1.918+1.2581 

1.891+1.2721 

1.725+1.2261 

0.675+0.17741 

LINEAR S'XMMETRICAL TORSION OF RECTANGULAR WING OF AR = 3 (a = O) 

k 

y 

c 

CL 


L 

boS 

2pU^oCt^ 

2 pU^o< 3 - 

0 

0 

0 

0 

0 

O.k 

0.8 

1.0 

1 

1 

1 

1 

00 

-0.352 

-0.596 

-1 

0 

1.257 

2.515 

5.142 

0 .k 55 

0.81k 

1.016 

0 

0.167 

0.167 

0.167 

0.167 

0 

o.k 

0.8 

1.0 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

oo -ooi 

-0.144+0.1201 

-0.365+0.1411 

-0,750+0.1051 

0 

0.971-0.0521 

1.942-0.1041 

2.427-0.1501 

0.406-0.0201 

0.777+0.0851 

1.001+0.1741 

0.044+0.0031 

0.555 

0.555 

0.555 

0.355 

0 

O.k 

0.0 

1.0 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

00 - 0^1 

-0.076+0.0821 

-0.262+0.0981 

-0.622+0.0101 

0 

0.555+0.1261 

1.707+0.2521 

2.153+0.5151 

0.540+0.0051 

0.741+0.2151 

1.007+0.5881 

0.175+0.0221 

0.667 

0.667 

0.667 

0.667 

0 

O.k 

0.8 

1.0 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.569-0.1301 

00 - ooi 

-0.044+0.0561 
-0.175+0.0121 
-0.459-0. 1851 

0 

0.770+0.4941 

1.559+0.9871 

1.924+1.2541 

0.255+0.0861 

0.699+0.5211 

1.095+0.8741 

0.675+0.1761 

LINEAE SYMMETRICAL TORSION OF RECTANGULAR WING OF AR = 6 (a = O) 

0 

0 

0 

0 

0 

O.k 

0.8 

1.0 

1 

1 

1 

1 

-0.2087 

-0.4290 

-1 

0 

1.257 

2.513 

5.142 

0.423 

0.994 

1.455 

0 

0.167 

0.167 

0.167 

0.167 

0 

O.k 

0.8 

1.0 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

00 - 00 i 

-0.049+0.0691 

-0.215+0.1281 

-0.750+0.1051 

0 

0.971-0.0521 

1.942-0.1041 

2.427-0.1501 

0.292-0.1021 

0.902+0.0501 

1.575+0.1731 

0.045+0.0031 

0.355 

0.555 

0.555 

0.555 

0 

O.k 

0.8 

1.0 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

00- cO i 

-0.053+0.0401 

-0.157+0.0951 

wO. 622+0. 0101 

0 

0.855+0.1261 

1.707+0.2521 

2.155+0.5151 

0.250-0.0821 

0.504+0.1691 

1.525+0.4291 

0.173+0.0221 

0.667 

0.667 

0.667 

0.667 

0 

O.k 

0.8 

1.0 

0.569-0.1501 

0.569-0.1501 

0.569-0.1301 

0.569-0.1501 

00 +00 1 

-0.012+0.0151 

-0.079+0.0511 

-0.459-0.1851 

0 

0.770+0.4941 

1.559+0.9871 

1.924+1.2541 

0.158-0.0071 

0 . 84 i+ 0 . 516 i 

1.301+0.9851 

0.675+0.1771 
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TABLE VIII (continued.) 


NACA TH No. 1195 


LINEAR SYMMETRICAL BENDING FOB RECTANGULAR WING OF AR = 3 


k 

y 

c 




L 


boS 

2 pV%oh /bo 

2pU^oH/bo 

0.167 

0.167 

0.167 

0.167 

0 

o.i^ 

0.8 

1.0 

0,757-0.1881 

0.757-0.1881 

0.757-0.1881 

0,757-0.1881 

00 — ooi 

-o.ilA+o.iaoi 
-0. 565+0. iJ+li 
-0.750+0.1051 

0 

0.022+0.1591 

0.0I4I4.+0.5171 

0.055+0,5971 

0.009+0.0671 
-0.005+0.1281 
-0.015+0.1651 
- 0.0002+0. od+i 

I 

• • • • 

0000 

0 

O.k 

0,8 

1.0 

0.650-0.1751 

0.650-0.1751 

0.650-0,1751 

0.650-0.1751 

00 - ooi 

-0.076+0.0821 

-0.262+0.0981 

-0.622+0.0101 

0 

0.005+0.2721 

0. 007 + 0 . 5 lAi 

0.008+016811 

0. 017+0. iio6i 
-0.051+0.2581 
-0.075+0.5251 
-0.005+0.0291 

0.667 

0.667 

0.667 

0.667 

0 

O.ii. 

0.8 

1.0 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

OQ -QOi 

-0.0I44+0.0561 
-0,175+0.0121 
-0.14.59-0. 1851 

0 

-0.170+0.1+771 

-0.5142+0.9541 

-0.1+26+1.1921 

-0.005+0.1581 
- 0.201+0.4401 
-0.561+0.6641 
-0.051+0.2501 

LINEAR SYMMETRICAL BENDING FOR RECTANGULAR WING OF AR = 6 

0.167 

0.167 

0.167 

0.167 

0 

o.k 

0.8 

1.0 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

0,757-0.1881 

00 - °oi 

-0.014.9+0.0691 

-0.215+0.1281 

-0.750+0.1051 

0 

0. 022+0. 1S91 
0.044+0.5171 

0.055+0.5971 

0.021+0,0471 

0.007+0.1481 

-0.010+0.2271 

-0.0002+0.004i 

0.555 

0.555 

0.555 

0.555 

0 

o.k 

0.8 

1.0 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

oO - 00 i 

-0.055+0.014.01 

-0.157+0.0951 

-0.622+0.0101 

0 

0.005+0.2721 

0.007+0.5441 

0.008+0.6811 

0.059+0,0701 
-0.015+0.2591 
-0.071+0.4501 
-0 .002+0.0291 

0.667 

0.667 

0.667 

0.667 

0 

O.ij- 

0.8 

1.0 

0.569-0.1501 

0.569-0.1^01 

0.569-0.1501 

0.569-0.1501 

00 - ooi 

-0.012+0.0151 
-0.079+0.0511 
-O.l4.59-O. 1851 

0 

-0.170+0.4771 

-0.541+0.9541 

-0.426+1.1921 

0.052+0.0821 

-0.185+0.4671 

-0.595+0.6221 

-0.042+0.2501 

PARABOLIC SYMMETRICAL BETOINa FOR RECTAWOULAR WING OF AR = 3 


k 

y 

C 


l (^) 


L 


"bo® 

2pU%oh^o 

2pU^oH^o 

0.167 

0.167 

0.167 

0.167 

0 

o.k 

0.8 

1.0 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

0.757-0.1881 

^ -001 
0.055+0.1061 
-0.592+0.114,51 
-0.750+0.1051 

0 

0.009+0.0651 

0.055+0.2541 

0.055+0.5971 

0.006+0.0561 
0.000+0.0661 
-0.015+0.1251 
-0.0002+0. 004 i 

0000 

• • • • 

J 

0 

Si 

1,0 

0.650-0.1751 

0.650-0.1751 

0.650-0.1751 

0,650-0.1751 

oO - ool 

0.079+0.0721 

-0.287+0.0971 

-0.622+0.0101 

0 

0.001+0.1091 

0.005+0.4561 

0.008+0.6811 

0.012+0.0591 

-0.011+0.1221 

-0.060+0.2451 

-0.005+0.0291 

0.667 

0.667 

0.667 

0.667 

0 

o.k 

0.8 

1.0 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

0.569-0.1501 

CP - toi 

0.069+0.0501 

-0.188+0.0061 

-0.14.59-0.185i 

0 

-0.068+0.1911 

-0.272+0.7651 

-0.426+1.1921 

0.005+0.0811 

-0.085+0.2141 

-0.280+0.5111 

-0.051+0.2501 
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PARABOLIC SYMMETRICAL BENDING FOR RECTANGULAR WING OF AR = 6 


0.167 

0 

0.757-0.1881 

00 - °ol 

0 

0.019+0.0251 

0.167 

O.k 

0.757-0.1881 

0.105+0.0201 

0.009+0.0651 

0.007+0.0721 

0.167 

0.8 

0.757-0.1881 

-0.259+0.1531 

0.055+0.2541 

-0.010+0.1741 

0.167 

1.0 

0.757-0.1881 

-0.750+0.1051 

0,055+0.3971 

-0.0002+0.0041 

0.355 

0 

0.650-0.1751 

• 00 - ooi 

0 

0.025+0.0551 

0.555 

O.k 

0.650-0.1751 

0.091-0.0051 

0.001+0.1091 

0.002+0.1241 

0.555 

0.8 

0.650-0.1751 

-0.155+0.0951 

0.005+0.4561 

-0.059+0.3531 

0.555 

1.0 

0.650-0.1751 

-0.622+0.0101 

0.005+0.6511 

-0.0002+0.0291 

0.667 

0 

0.569-0.1501 

00 +^i' 

0 

0.020+0.0361 

0.667 

o.k 

0.569-0.1301 

0.057+0. 0011 

-0.068+0.1911 

-0.068+0.2101 

0.667 

0.8 

0.569-0.1301 

-0.086+0.0271 

-0.272+0.7651 

-0.509+0.6481 

0.667 

1.0 

0.569-0.1301 

-0.k59-0.1851 

-0.426+1.1921 

-0.042+0.2301 

LINEAR AHTISYMMETRICAL TORSION OF BECTAIWULAR WING AR = 3 


7 




L 

k 


c 


SpU^oJ 

2 pU^oI 

0 

0.4 

1 

-0.490 

1.257 

0.641 

0 

0.6 

1 

-0.52jB 

1.885 

0.890 

0 

0.8 

1 

-0.584 

2.515 

1.947 

0 

1.0 

1 

-1 

5.142 

0 

0.167 

0.4 

0.757-0.1881 

-0.253+0.1651 

0.971-0.0521 

0.656+0.1291 

0.167 

0.6 

0.757-0.1881 

-0.289+0.1621 

1.456-0.0781 

0.887+0.1821 

0.167 

0.8 

0.757-0.1881 

-0.342+0.1591 

1. 942-0. loki 

1.049+0.2241 

0.167 

1.0 

0.757-0.1881 

-0.750+0.1051 

2.427-0.1501 

0 .o 44 '*'^ .0031 

0.555 

o.k 

0.650-0.1751 

-0.149+0.1501 

0.855+0.1261 

0.659+0.2591 

0.555 

0.6 

0.650-0.1751 

-0.185+0.1281 

1.260+0.1891 

o.895'*"0.574i 

0.555 

0.8 

0.650-0.1751 

-0.255+0.1221 

1.707+0.2521 

1.069+0.4611 

0.555 

1.0 

0.650-0.1751 

- 0 . 622 + 0.0101 

2.135+0.5151 

0.175+0.0221 

0.667 

o.k 

0.569-0.1501 

-0.072+0.0611 

0.770+0.4941 

0.645+0.5371 

0.667 

0.6 

0.569-0.1501 

-0.108+0.0541 

1.155+0.7411 

0.917+0 ..7741 

0.667 

0.8 

0.569-0.1501 

-0.l46+0.0k0i 

1.559+0.9671 

1.158+0.9661 

0.667 

1.0 

0.569-0.1301 

-0.459-0.1851 

1.924+1.2541 

0.675+0.1781 


RIGID ROLLING OF RECTANGULAR WING OF AR = 3 

0.167 

O.k 

0.757-0.1881 

-0.255+0.1651 

0.022+0.1591 

-0.013+0.1061 

0.167 

0.6 

0.757-0.1881 

-0.269+0.1621 

0.055+0.2561 

-0.018+0.1471 

0.167 

0.8 

0.757-0.1881 

-0.542+0.1591 

0.044+0.5171 

-0.023+0.1741 

0.167 

1.0 

0.757-0.1881 

-0.750+0.1051 

0.055+0.5971 

-0.0002+0 ,oo4i 

0.555 

O.k 

0.650-0.1751 

-0.149+0.1501 

0.005+0.2721 

-0.051+0.2101 

0.555 

0.6 

0.650-0.1751 

-0.165+0.1281 

0.005+0.4081 

-0.076+0.2951 

0.555 

0.8 

0.650-0.1751 

-0.253+0.1221 

0.007+0.5441 

-0.094+0.3511 

0.555 

1.0 

0.650-0.1751 

-0.622+0.0101 

0.006+0.6811 

-0.005+0.0291 

0.667 

O.k 

0.569-0.1501 

-0.078+0.0611 

-0.170+0.4771 

- 0 . 222 + 0.4111 

0.667 

0.6 

0.569-0.1501 

-0.108+0.0541 

-0.255+0.7151 

-0.523+0.5791 

0.667 

0.8 

0.569-0.1501 

-0.146+0.0401 

-0.541+0.8541 

-0.408+0.7081 

0,667 

1.0 

0.569-0.1501 

-0.459-0. 1651 

-0.426+1.1921 

-0.042+0.2301 

1.555 

o.k 

0.526-0.0811 

-0.054-0.0011 

-0.981+0.8811 

-0.979+0.8241 

1.555 

0.6 

0.526-0.0811 

-0.042-0. 0151 

-1.471+1.5221 

-1.454+1.2151 

1.555 

0.8 

0.526-0.0811 

-0.055-0.0521 

-1.961+1.7651 

-1.853+1.5861 

1.555 

1.0 

0.526-0.0811 

-0.090+0.4141 

-2.452+2.2051 

-0.716+1.8251 
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LINEAR SYMMETRICAL TORSION OF ELLIPTICAL WING OF AR = 3 fa = O) 


k 

.7 

^ 0 ® 

c 


E(-) 

2 pU^lDQa 

L 

2 pU^t>oa 

0 

0 

0 

0 

O.k 

0.8 

1 

1 

1 

oo 

-0.1419 

-0.555 

0 

1.152 

1.508 

0.591 

0.669 

0.702 

0.212 

O.I9I4- 

0.127 

0 

O.k 

0.8 

0.719-0.1881 

0.755-0.1891 

0.799-0.1811 

00 - ooi 

-0. 181+0. 11I4.1 
-0.514.7+0.114.11 

0 

0.865-0.0251 

1.222-0.1021 

O.5I49-O.OI2I 

0 . bkk+O . 0881 
0.685+0.0781 

O.I4.2L. 

0.588 

0 . 25 k 

0 

O.k 

0.8 

0.618-0.1621 

0.629-0.1671 

0.691-0.1851 

00 - ooi 

0. 029+0. 05ki 
-0.2714+0.12141 

0 

0.762+0.1721 
1. 077+0. 0I461 

O.29I4+O .0571 
0.650+0.2501 
0.642+0.1851 

0 . 8[(.7 

0.776 

0.508 

0 

O.k 

0.8 

0.551-0.1151 

0.556-0.1191 

0.596-0.1501 

00 - OOl 

-0. 060-0. 0111 
-0.177+0.0661 

0 

0.6911+0.5591 

0.956+0.5861 

0. 175+0. llili 

0 . 650+0 . R20i 

0 .66I4.+O .4l8i 

L 

INEAB SYMMETRICAL TORSION OF ELLIPTICAL WING OF AR = 6 ( 

a = 0) 

0 

0 

0 

0 

1 

1 

1 

+ CO 

-0.26k 

-0.576 

0 

1.152 

1.508 

0.588 

0.849 

o. 9 ki 

0.212 

0 . 19 ^ 4 - 

0.127 

0 

O.k 

0.8 

0.719-0.1881 

0.755-0-1891 

0.799-0.1811 

00 - ooi 

-0.080+0.0761 
-0. 209+0. 12I4I 

0 

0.865-0.0251 

1.222-0.1021 

0.294-0.0711 
0 .764+0 .0561 
0.895+0.0651 

0.14.2k 

0.588 

0.2514 

0 

O.k 

0.8 

0.618-0.1621 

0.629-0.1671 

0.691-0.185 

00 - ooi 

-0 .0l|.6+0 .0[|.5i 
-0.114.5*^0.1121 

0 

0.762+0.1721 
1. 077+0. 0I461 

0.212-0.0441 

0.700+0.2151 

0.857+0.1871 

o. 8 k 7 

0.776 

0.508 

0 

O.k 

0.8 

0.551-0.1151 

0.556-0.1191 

0.596-0.1501 

00 - ooi 

-0. 026-0. ooki 
-0. 088+0. 06I4I 

0 

0.69I4+0.5591 

0.956+0.5561 

0.100+0.0251 

0.680+0.51441 

0.800+0.4411 


LINEAR SYMMETRICAL BENDING OF ELLIPTICAL WING AR = 3 \ 

k 

7 

c 


l (^) 

2pU%oiiAo 

L 

2pI^^)oh^o 

0.212 

0.1914- 

0.127 

0 

o.k 

0.8 

0.719-0.1881 

0.755-0.1891 

0.799-0.1811 

- 1 

-0.180+0.1191 

-0.5147+0.1591 

0 

0.022+0.1791 

0.058+0.2551 

0.006+0.0741 

-0.007+0.1551 

-0.007+0.11441 

0.14.2k 

0.588 

0.25I4. 

0 

o.k 

0.8 

0.615-0.1621 

0.629-0.1671 

0.691-0.1851 

- 1 

-0.115+0.0751 

-0.277+0.1161 

0 

-0.015+0.^071 

0.057+0.141421 

0.016+0.0701 

-0.050+0.2511 

-0.057+0.2641 

0.8I47 

0.776 

0.508 

0 

0.14. 

0.8 

0.551-0.1151 

0.556-0.1191 

0.596-0.1501 

- 1 

-0.057+0. 0011 
-0.177+0.0651 

0 

-0. 265+0. 5 k 5 i 
-0.15L4+0.7621 

-0.059+0.1651 

-0.279+0.4971 

-0.2114+0.5561 


NATIONAL ADVISORY 


COMMITTEE FOR AERONAUTICS. 



NACA IN No. 1195 


TABLE VIII (continued) 


81 


LINEAR SYMMETRICAL BENDING OF ELLIPTICAL WING OF AR = 6 


0.212 

0 . 19 ii 

0.127 

0 

0.4 

0.8 

0.719-0.1881 

0.755-0.1891 

0.799-0.1811 

- 1 

-0.078+0.0801 

-0.211+0.1251 

0 

0.022+0.1791 

0.058+0.2551 

0.021+0.0611 

0.005+0.1601 

-0.002+0.1881 

0 .k 2 k 

0.588 

0.25k 

0 

0.4 

0.6 

0.618-0.1621 

0.629-0.1671 

0.691-0.1851 

- 1 

-0.042+0.0411 

-o.i46+o.iioi 

0 

-0.015+0.5071 

0.057+0.4421 

0.054+0.0851 

-0.054+0.2871 

-0.054+0.5491 

0.847 

0.776 

0.508 

0 

0.4 

0.8 

0.551-0.1151 

0.556-0.1191 

0.596-0.1501 

- 1 

-0.021+0.0021 

-0.087+0.0651 

0 

-0.265+0,5451 

-0.154+0.7621 

0.011+0.0961 

-0.266+0.5241 

-0.216+0.6521 

PARABOLIC SYMMETRICAL BENDING OF ELLIPTICAL WING OF AE = 3 

k 

7 

c 

0 

L(2) 

L 

tos 

2pU^oh^o 

2pU^oh/6o 

0.212 

0.194 

0.127 

0 

0.4 

0.8 

0.718-0.1881 

0.755-0.1891 

0.799-0.1811 

oO — oOl 

-0.050+0.1051 

-0.572+0.1441 

0 

0.009+0.0711 

0.050+0.2041 

0.005+0.0591 

-0.001+0.0691 

-0.007+0.1091 

0.424 

0.588 

0.254 

0 

0.4 

0.8 

0.618-0.1621 

0.629-0.1671 

0.691-0.1851 

oO - ^ 1 

0.029+0.0541 

-0.274+0.1241 

0 

-0.005+0.1251 

0.029+0.5551 

0.005+0.0651 

-0.016+0.1281 

-0.054+0.2151 

0.847 

0.776 

0.508 

0 

0.4 

0.8 

0.551-0.1151 

0.556-0.1191 

0.596-0.1501 

oO - oO i 

0.059+0.0121 

-0.187+0.0611 

0 

-0.105+0.2171 

-0.107+0,6101 

-0.010+0.0851 

-0.116+0.2551 

-0.169+0.4181 

PARABOLIC SYMMETRICAL BENDING OF ELLIPTICAL WING OF AR = 6 

0.212 

0.194 

0.127 

0 

0.4 

0.8 

0.719-0.1881 

0.735-0.1891 

0.799-0.1811 

OO - oO i 

0.049+0.0581 

-0.228+0.1281 

0 

0.009+0.0721 

0.050+0.2041 

0.012+0.0501 

0.005+0.0761 

-0.005+0.1461 

0.424 

0.588 

0.254 

0 

0.4 

0.8 

0.618-0.1621 

0.629-0.1671 

0.691-0.1851 

CiO - CO 1 

0.049+0.0051 

-0.157+0.1151 

0 

-0.005+0.1251 

0.029+0.5551 

0.021+0.0591 

-0.006+0.1551 

-0.028+0.2751 

0.847 

0.776 

0.508 

0 

0.4 

0.8 

0.551-0.1151 

0.556-0.1191 

0.596-0.1501 

CJO - 00 1 
0.038-0.0151 
-0.087+0.0591 

0 

-0.105+0.2171 

-0.107+0.6101 

0.012+0.0421 

-0.101+0.2521 

-0.145+0.5211 
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LINEAR ANTISIMMETRICAL TORSION OF ELLIPTICAL WING OF AR = 3 


k 

7 

toS 

c 


l(=) 

2pU^oJ 

L 

2pU%oj 


1^0 = ° . 


0 

0 

0 

O.k 

0.6 

0.8 

1 

1 

1 

-0.571 

-0.571 

-0.571 

1.152 

1.508 

1.508 

0.494 

0.647 

0.647 


k = 0.212 
0 


0 . 19 k 

0.170 

0.127 

o.k 

0.6 

0.8 

0.755-0.1891 

0.754-0.1881 

0.799-0.1811 

-0.289+0.1521 
-0 . 504+0 .1551 

-0.505+0.1621 

0.865-0.0251 

1.162-0.0591 

1.222-0.1011 

0.515+0.1191 

0.685+0.1561 

0.747+0.1141 



k^ = 0 . k 2 h 


0.588 

0.55? 

0.25!+ 

o.k 

0.6 

0.8 

0.629-0.1671 

o.6i^.7-o.l7kl 

0.691-0.1851 

-0.185+0.0961 

-0.197+0.1141 

-0.200+0.1501 

0.762+0.1721 
1 . 021 + 0 . 159 i 
1.077+0.0461 

0.550+0.2421 

0.695+0.2801 

0.747+0.2551 



k^ = o. 8 k 7 


0.777 

0.678 

0.509 

o.k 

0.6 

0.8 

0.562-0.1241 

0.569-0.1501 

0.569-0.1501 

-0. 099+0. OOlil 
-0.114+0.0501 
-0.121+0.0951 

0.705+0.5561 

0.924+0.6071 

0.957+0.5861 

0.587+0.5171 

0.756+0.5951 

0.758+0.4801 


RIGID ROLLING OF ELLIPTICAL WING OF AR = 3 


k„ = 0.212 

0 


0 . 19 i 4 - 

0.170 

0.127 

O.k 

0.6 

0.8 

0.755-0.1891 

0.754-0.1881 

0.799-0.1811 

-0.289+0.1561 

-o.5oL(.+o.158i 

-0.504+0.1621 

0.022+0.1791 

0.055+0.2411 

0.058+0.2551 

-0.016+0.1081 

-0.018+0.1441 

-0.014+0.1581 



k^ = 0.l).2l+ 


0.588 

0.35? 

0.25^^ 

O.k 

0.6 

0.8 

0.629-0.1671 

0 . 6 l 4 - 7 - 0 . 17 ki 

0.691-0.1851 

-0. 182+0. lOkl 
-0.196+0.1161 
-0. 200+0. ik8i 

-0.015+0.5071 

0.005+0.4141 

0.057+0.4411 

-0.064+0.2181 

-0.075+0.2881 

-0.058+0.5151 


k^ = o. 8 k 7 


0.777 

0.678 

0.507 

o.k 

0.6 

0.8 

0.562-0.1241 

0.569-0.1501 

0.596-0.1501 

-0.098+0.0161 

-0.115+0.0571 

-0.120+0.0891 

-0.258+0.5491 

-0.268+0.7271 

-0.154+0.7621 

-0.274+0.4541 

-0.515+0.5631 

-0.246+0.6081 



= 1.695 


1.555 

1.556 
1.017 

O.lf 

0.6 

0.8 

0.520-0.0721 

0.525-0.0801 

0.559-0.0991 

-0.005-0.0571 
-0.050-0. 0551 
-0.056+0.0261 

-1.576+1.0151 
-1.527+1.5451 
-1.046+1 . 5771 

-1.304+1.0061 

-1.438+1.2651 

-1.113+1.2551^ 
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UNIFOBM PITCHIHG OF ELLIPTICAL WHG OF AR = 3 


k 

7 

c 



L 


2pU®b^J 

2pU®b^J 

0 

0 

0 

0 

1 

1 

1 

-o.koo 

-o.koo 

-o.koo 

2.878 

1.727 

1.885 

1.727 

1.056 

0.212 

0 . 19 U 

0.127 

0 

O.k 

0.8 

0.719-0.1881 

0.733-0.1891 

0.799-0.1811 

-o.ik8+o.ii8i 

-0.162+0.1221 

-0.227+0.1511 

2.522+0.0061 

2.161-0.0591 

1.527-0.1261 

1.816+0.5271 
1. 659+0. 2k81 
1.084+0.0951 

O.U2k 

0.388 

0.25k 

0 

0.618-0.1621 

0.629-0.1671 

0.691-0.1851 

-0.078+0.0631 

-0.085+0.0811 

-0.151+0.1181 

2 .ok 9 +o. 569 i 

i. 9 ok+o.k 30 i 

i. 3 k 7 +o. 057 i 

1.760+0.7161 
1.620+0.6181 
1 .03k+o.2kki 

o. 8 k 7 

0.776 

0.508 

0 

o.k 

0.8 

0.551-0.1131 

0.557-0.1191 

0.596-0.1501 

-0.057+0.0011 

-o.ok9+o.oi3i 

-0.101+0.0791 

1.880+1.7101 

1 . 755 + 1-5571 

1. 195+0. k85i 

1.761+1.6651 

1.579+1.5801 

0.967+0.5851 

1.695 

1.555 

1.017 

0 

o.k 

0.8 

0.518-0.0671 

0.520-0.0721 

0.559-0.0991 

0.005-0.0181 

-0.006-0.0221 

-0.081+0.0201 

1.861+3.8291 

1.657+5.1891 

1.110+1.2881 

1.925+5.7871 
.1.687+5.1151 
0.959+1 - 2 k 7 i 

UNIFORM TRANSLATION OF ELLIPTICAL WING AR = 3 

0.212 

O.I9U 

0.127 

0 

o.k 

0.8 

0.719-0.1881 

0.755-0.1891 

0.799-0.1811 

-0.lk8+0.1211 
-0. 162+0. 12kl 
-0.228+0.1251 

0. 055+0. k78i 
0. 056+0. kk7i 
o.ok 7 +o. 5 i 9 i 

-0.026+0.5801 

-0.020+0.5481 

-0.005+0.2281 

o.k2k 

0.388 

0.25k 

0 

o.k 

0.8 

0.618-0.1621 

0.629-0.1671 

0.691-0.1851 

-0.076+0.0701 

-0.082+0.0561 

- 0 .i 5 k+o.i 07 i 

-0.067+0.8251 

-0.055+0.7671 

o.ok 6 +o. 552 i 

-0.160+0.7211 

-0.158+0.6671 

-0.040+0.4291 

0 . 8 k 7 

0.776 

0.508 

0 

Si 

0.551-0.1151 

0.557-0.1191 

0.596-0.1501 

-0.056+0.0081 

-o.ok 7 +o.oi 7 i 

-0.105+0.0651 

-0. 828+1. k66i 
-0.657+1.5581 
-0.167+0.9521 

-0.849+1.5691 

-0.698+1.2451 

-0.271+0.7671 

1.695 

1.555 

1.017 

0 

O.k 

0.8 

0.518-0.0671 

0.520-0.0721 

0.559-0.0991 

0.005-0.0121 

-0.006-0.0161 

-0.075+0.0081 

-k- 154+2.7571 

- 3 .k 56 + 2 . 557 l 

-1.508+1.7211 

-k. 088+2. 7731 

-5.557+2.5061 

- 1 . 355 + 1.4811 


LINEAR ANTISYMMETRICAL TORSION OF AN ELLIPTICAL WING OF AR = 3 
FOR HIGHER VALUES OF kg 

V 1.695 

1.555 

1.556 
1.017 

O.k 

0.6 

0.8 

0.520-0.0721 

0.525-0.0801 

0.559-0.0991 

-0. 001-0. ok9i 
-0. 030-0. ok5i 

-0.075+0.0161 

o. 665 +r. 277 i 

0. 874+1. k5 81 
0.888+1.0501 

0.706+1.2201 

0.676+1.5401 

0.765+0.9981 


ko= 2 . 5 k 2 


2.550 

2.03k 

1.525 

Si 

0.8 

0.511-0.0521 

0.515-0.0571 

0.521-0.0751 

0.059-0.0001 

0.031-0.0311 

-0.057-0.0301 

0 . 6 k 8 + 1 . 967 l 
0. 861+2. 25ki 
0.869+1.6591 

0.695+2.0201 

0.954+2.2541 

0 . 847 + 1.5521 
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UMIFOEM SPAIWISE AILEBON DEFLECTION OF A FULL-SPAN AILERON 
OF ELLIPTICAL WIND OF AR = 3 (e = 0 , I = 0.2 'bA’o) 


k 

1 

c 

Op 


L 

boB 

2pU^oP 

2pU^oP 

0 

0 

1 

-0.1+00 

2.589 

1.652 

0 

0.1+ 

1 

-0.1+00 

2.375 

I. 5 li+ 

0 

0.8 

1 

-0.1+00 

1.555 

0:991 

0 .14.2I1. 

0 

0.618-0.1621 

-0.080+0.0621 

" 1 . 661 ++ 0.5791 

1.1+11+0.1+811 

0.589 

0.4 

0.629-0.1671 

-0. 090+0. 071+1 

1.551++0.2821 

1.29IL+0.I+001 

0.25^ 

0.8 

0.691-0.1851 

-0.150+0.1211 

1.111+0.0201 

0.81+6+0. 1671 

0 . 81+7 

0 

0.551-0.1151 

-O.Ol+l+O.OOll 

1.569+1.2001 

1.26k+l.ll+ll 

0.776 

0.1+ 

0.556-0.1191 

-0.050+0.0091 

1.317+0.971+1 

1 .186+0 .9521 

0.508 

0.5 

0.596-0.1501 

-0.098+0.0851 

0.958+0.5291 

0.761+0.4061 


UNIFORM SPANWISE AILERON DEFLECTION 

OF A FULL-SPAN AILERON 



OF ELLIPTICAL WING OF AE = 3 (e 

= 0, 1 = 0.2 b/bjj) 


y 



Mp(^)(e) 

Mp (e) 

k 

boS 

C 

Op 

2pU^q2^ 


0 

0 

1 

-0.1+00 

0.187 

0.124 

0 

o.k 

1 

-0.1+00 

0.157 

O.IOI+ 

0 

0.8 

1 

-0.1+00 

0.067 

O.OI+5 

O.Il2L 

0 

0.618-0.1621 

-0.080+0.0621 

0.097+0,1911 

0.081+0.1971 

0.589 

o.k 

0.629-0.1671 

-0. 090+0. 07I+1 

0.087+0.1551 

0.075+0.1611 

0.25L1. 

0.8 

0.691-0.1851 

-0.150+0.1211 

O.OI+7+O.O57I 

0.057+0.01+21 

0 . 81+7 

0 

0.551-0.1131 

-O.Ol+l+O.OOll 

-0.025+0.1+081 

-0.029+0.1+061 

0.776 

o.k 

0.556-0.1191 

-0.050+0.0091 

0.005+0.5121 

-0 . 005+0 .5101 

0.508 

0.8 

0.596-0.1501 

-0.098+0.0831 

0. 025+0. 08I+1 

0.021+0.0871 


UNIFORM AILERON HINGE DEFLECTION OF A HALF-SPAN AILERON OF AN ELLIPTICAL WING 


OF AR = 6 


(• ' 5 - 


I = 0.1 


-) 

bo/ 
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0 

0 

0 

0 

0 


0.212 

0 .i 9 i; 

0.170 

0.127 

0.092 


0 .1+214. 
0.589 
0.55? 
0.251+ 
0.185 


bns 


n 

0.8 

0.9 


0 

o.k 

0.6 

0.8 

0.9 


O.li 

0.6 

0.8 

0.9 


0.719-0.1881 
0.735-0.1891 
0.75f+-0.i88i 
0.799-0.1811 
0.81+5-0. 1681 


0.618-0.1621 
0.629-0.1671 
0. 61+8-0. 17I+1 
0.691-0.1851 
0.71+1-0.1881 


00 


- 0 . 1 + 91 + 

-0.560 

-0.259 


00 - ool 
00 — oO 1 

-0.275+0.1261 

-0.188+0.1501 

-0.116+0.0911 


CM 3 - ool 

00 - 

-0.187+0.0871 

■0.107+0.1271 

-O.OI+I1.+O.II+81 


_L^ 

2pU^nP 


0 

0 

1.551^ 

0 ^ 7 ^ 


0 

0 

1.187-0.1821 

0.950-0.11+81 

0.629-0.1021 


0 

0 

1.052-0.0701 

0.813-0.1001 

0.557-0.0991 


gpU^nP 


0.080 

0.255 

0.786 

0.758 

0.566 


0.056-0.0561 

0.175-0.1211 

0.752-0.0101 

0.707-0.0071 

0.51+1-0.0561- 


0.006-0.0571 

0.108-0.0261 

0.726+0.0511 

0.677+0.0561 

0.519+0.0101 
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UNIFOEM AILERON KENGE DELFECTION OF A HALF-SPAN AILERON OF AN ELLIPTICAL WING 


OF AR = 6 (^e = -, I = 0.1 —') 

\ 2 


k 

y 

c 

03 

Mp^^)(e) 

M 3 (e) 
2pU®bo®3 



ko= 0 



0 

0 

0 

0.6 

0.8 

0.9 

1 

1 

1 

-0.494 

-0.560 

-0.239 

0.050 

0.015 

0.005 

0.025 

0.015 

0.005 

k =0.212 
o 

0.170 

0.127 

0.092 

0.6 

0.8 

0.9 

0.754*0.1881 

0.799-0.1811 

0.845-0.1681 

-0.275+0.1261 

-0.188+0.1501 

-0.116+0.0911 

0.026+0.0051 

0.014+0.0021 

0.005+0.0001 

0.022+0.0071 

0.012+0.0051 

0.005+0.0011 

kQ= 0.1+214. 

0.559 

0.254 

0.185 

0.6 

0.8 

0.9 

0 . 648 - 0.1741 

0.691-0.1851 

0.741-0.1881 

-0.187+0.0871 

-0.107+0.1271 

-0.044+0 *1401 

0.024+0.0131 

0.013+0.0051 

0.005+0.0011 

0,021+0.0141 

0.012+0.0061 

0.005+0.0011 

ROLLING-MOMENT COEFFICIENTS Cp FOR A RECTANGULAR WING OF AR » 3 

^0 

Due to rigid rolling 

Due to linear antleymmetrlced 
torsion 

Or/h 


Cp/a 

Cr(*)/5 

0 

0.167 

0.335 

0.667 

1.353 

0 

-0.008+0.0621 

-0.053+0.1261 

-0.142+0.2621 

-0.676+0.6561 

0 

0.018+0.1521 

0.005+0.2271 

-0.142+0.3971 

-0.817+0.7541 

0.372 

0.575+0.0781 

0.588+0.1611 

0.459+0.5451 

0.620+0.8551 

1.047 

0.809-0.0441 

0.711+0.1051 

0.641+0.4111 

0.608+0.9801 

R( 

DLLING-+10MENT COEFFICIENTS Cp FOR AN ELLIPTICAL WING OF AR = 3 

ko 

Due to rigid rolling 

Due to linear antisynmietrical 
torsion 

«CR/ 4 h 




0 

0.212 

O.L2I4. 

0.814-7 

1.695 

2.542 

0 

-0.006+0.0581 

-0.026+0.1161 

-0.111+0.2521 

-0.511+0.14.861 

-1.251+0.7551 

0 

0.015+0.0951 
0.00 6+0. i64i 

-0.082+0.2871 

-0.515+0.5241 

-1.256+0.7641 

0.250 

0.276+0.0491 

0.279+0.1091 

0.288+0.2111 

0.519+0.4691 

0.541+0.7561 

0.502 

0.458-0.0291 

0,404+0.0441 

0.565+0.1991 

0.540+0.4911 

0.555+0.7681 
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LIFT COEFFICIENTS Cl FOB AN ELLIPTICAL WING OF AE = 3 


ko 

Due to uniform translation 

Due to voilform pitching 

jiCt j — 
^ 8h 

/g- 

GO 

P 1 

0 

0.212 

0.k2l|. 

0.847 

1.695 

0 

-0. 001+0. 29ki 
-0.10^+0.5661 
-0.554+1.0561 
-2.657+2.0801 

0 

0.048+0.5891 

-0.015+0.6721 

-0.494+1.1721 

-2.665+2.1651 

1.4B0 

1.420+0.2121 

1.556+0.4561 

1.557+1.0891 

1.585+2.4211 

2.467 

1.867-0.0691 

1.665+0.2971 

1.495+1.0751 

1.451+2.4681 


RECTANGULAE V/ING OF AR = 13.5 WITH INDICATED MODES (k^ = 0 . 097 ) 


Mode 

y/boS 

Torsion 

First 

bending 

Second 

tending 

0 

0 

0 

0 

0.4 

0.586 

0.250 

0.685 

0.6 

0.809 

0.461 

0.589 

0.8 

0.951 

0.725 

-0.068 

1.0 

1.0 

1.0 

-1.0 


Torsion mode 

k 

y/bo^ 

c 


E(^) 

L 

2 pU%nI 

0 0 0 0 0 
00000 

0 

0.4 

0.6 

0.8 

1.0 

0.857-0.1701 

0.857-0.1701 

0.857-0.1701 

0.857-0.1701 

0.857-0.1701 

C 3 Q - 1 

-0.045+0.0481 

-0.069+0.0711 

-0.109+0.0891 

-0.855+0.1221 

0 

1.565-0.1511 

2.159-0.2091 

2-558-0.2451 

2.669-0.2581 

0.502-0.1401 

1.472-0.0701 

1.968-0.0451 

2.190-0.0081 

0.015-0.1021 

FIRST BENDING 

0.097 

0.097 

0.097 

0.097 

0.097 

0 

0.4 

0.6 

0.8 

1.0 

0.857-0.1701 

0.857-0.1701 

0.657-0.1701 

0.857-0.1701 

0.857-0.1701 

OO - <501 

0.056+0.0011 

-0.046+0.0641 

-0. 124+0. 096I 
-0.855+0.1221 

0 

+0.009+0.0581 
+0.017+0.1171 
0. 027+0. i 84 i 
0.0^7+0.2541 

0.009+0.0151 
0.008+0.0621 
0.008+0. 1111 
0.006+0.1571 
+0.000+0.0011 

SECOND BENDING 

0.097 

0.097 

0.097 

0.097 

0.097 

0 

0.4 

0.6 

0.8 

1.0 

0.857-0.1701 

0.857-0.1701 

0.857-0.1701 

0.857-0.1701 

0.857-0.1701 

- 1 

-0.155+0. 1011 
-0.176+0.1051 
-0.619+0.5091 
-0.855+0.1221 

0 

0.025+0.1751 
0.022+0.1501 
-0. 005-0.0171 
-0.057-0.2541 

.0.011+0.0261 

0.004+0.1421 

0.005+0.1181 

o.oo4-o.oo51 

0.000-0.0011 
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TABLE IX.- SPAN COPKECTIONS AT MIDSPAN FOR RIGID ELLIPTICAL WING 
IN TRANSLATION AND PITCHING. (ONE-POINT APPROXIMATION) 


ko 

kos 

C 

ah = Oct , = 0 

C + a 

AR 

0 

0 

1 + 01 

- 0.571 

0.429 



0,424 

0,5 

0.618 - 0.1621 

- 0.196 + 0.0921 

0.422 r .. 0.0701 



0.847 

1.0 

0.551 - 0.1131 

- 0.114 + 0.0081 

0.437 - 0.1051 

\ 

> 1.5 

1,695 

2.0 

0.518 - 0.0661 

- 0.018 - 0.0441 

0.500 - 0.1101 

( 

3.390 

4.0 

0.505 - 0.0361 

+ 0.005 + 0.0071 

0.511 - 0.0291 

J 


0 

0 

1 + 01 

- 0.500 

0.500 



0.318 

0.5 

0‘.657 - 0.1771 

- 0.176 + 0.1161 

0.481 - 0.0611 



0.637 

1.0 

0.573 - 0.1341 

- 0.103 + 0.0481 

0.470 - 0.0861 


? 2 

1.273 

2.0 

0.528 - 0.0841 

- 0.042 - 0.0181 

0.486 - 0.1021 

1 

r 

2.546 

4.0 

0.509 - 0.0481 

+ 0.004 + 0.0171 

0.513 - 0.0311 

J 

1 

0 

0 

1 + 01 

- 0.400 

0,600 



0.212 

0.5 

0.719 - 0.1881 

- 0.154 + 0.1251 

0.565 - 0.0631 



0.424 

1.0 

0.618 - 0.1621 

- 0.079 + 0.0751 

0.539 - 0.0871 



0.847 

2.0 

0.551 - 0.1131 

- 0.043 + 0.0181 

0,507 - 0.0951 


1.695 

4.0 

0.518 - 0.0671 

- 0.010 - 0.0131 

0.508 - 0.0801 



2.542 

6.0 

0.509 - 0.0471 

+ 0.007 - 0.0051 

0.516 - 0.0521 

J 

I 

0 

0 

1-01 

- 0.250 

0.750 


1 

0.106 

0.5 

0.824 - 0.1751 

- 0.108 + 0.0961 

0.716 - 0.0781 



0.212 

1.0 

0.719 - 0.1891 

- 0.052 + 0.0771 

0.669 - 0.1111 


1 (■ 

0.424 

2.0 

0.618 - 0.1621 

- 0.021 + 0.0371 

0.597 - 0.1241 


( 

0.847 

4.0 

0.551 - 0.1131 

- 0.013 + 0.0091 

0.538 - 0.1031 



1.273 

6.0 

0.528 - 0.0841 

- 0.008 - 0.0001 

0.520 - 0.0841 


r 
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NACA TN Ho. 1195 


TABLE X.- IKIERMIllAHT TERMS AS A FUHCTIOH OF k FOR EXAMPIF I 







Determinant terms 
(tvo-dim.) 

/ 

Term 

Structural part of tern 

Aero, part of term 
(two-dlm.) 

\g * Sn * 8a 
and 






=n) 

k 0.333 

A 

1.3331 


C ? T * V 1 

-0.0096 - 0.77991 

1.3235 - 0.77991 - 0.25601 

B 

0.315U 



-3.0592 - 0.74371 

-2.7438 - 0 . 7437 " 

D 

0.3151* 



0.1274 + 0.19501 

0.4U28 + 0.19501 

£ 

0.8789 



0.9075 - 0.80171 

1.7864 - 0.8017 - 0.87890 

k » O.U 

A 

1.3331 



0.0350 - 0.62501 

1.3681 - 0.62501 - O.256CD 

B 

O . Zl^h 



-2.0U3I - 0.73441 

-1.7277 - 0.73441 

D 

O.315U 



0.1163 + 0.15631 

0.4317 + 0.15631 

£ 

0.8789 



0.6365 - 0.63751 

1.5154 - 0.63751 - 0.87890 

k » 0.5 

A 

1.3331 


»■!’“(?)’ <*• Vl 

0.0794 - 0.47831 

1.4125 - 0.47831 - O.256OR 

B 

0 . 315 i^ 



- 1.2414 - 0.67691 

-0.9260 - 0.67691 

D 

0.3154 



0.1051 + 0.11961 

0.4205 + 0.11961 

1 

0.8789 


(?)■»•■«.>! 

0.4227 - 0.48621 

1.3016 - 0.48621 - 0.87890 

k = 0.6 

A 

1.3331 


/ <u 1 

0.1920 [^) (1 . Ig^) 1 

0.1081 - 0.38591 

1.4412 - 0.38591 - 0.25600 

B 

0.3154 



-0.8208 - 0.61121 

-0.5054 - 0.61121 

D 

0.3154 



0.0980 + 0.09651 

0.4134 + 0.0961 

E 

0.8789 



0.3535 - 0.49271 

1.2513 - 0.4927 - 0.87890 
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TABLE XI 


SOLUTION OF SIMULTANEOUS EQUATIONS BY USE OF THE CEOUT METHOD OF REFERENCE 8 (EXAMPLE I) 







Llneeir torsion 

Parabolic bending 







Check 


Check 

Given Matrix 












1.0169 

-0.3837 

0.2702 

-0.2262 



0 

0.6772 



0 

0.6772 

- 0.05241 

0.07391 

-0.07851 

0.07901 




0 .02201 





.9401 

-.1323 

-.1299 

.2289 



0.4 

1.3068 



0.16 

1.0668 

-.05571 

.02621 

.04051 

-.08261 




-.07161 




-.05071 

.6672 

.4213 

-.0234 

-.3057 



.8 

1.5594 



.64 

1.3994 

-.06841 

-.12381 

.02451 

.11701 




-.05071 




-.05071 

.2966 

.4701 

. 723 ** 

.9896 



1.0 

3.4797 



1.0 

3.47971 

.09721 

-.05051 

-.18791 

-.30721 




-.44841 




.44841 

Supplementary Matrix 












0.9808 

-0.3837 

0.2702 

-0.2262 



0 

0.6772 



0 

0.6772 

.05051 

.07391 

-.07851 

.07901 




.02201 




.02201 

4.3261 


-.3793 

.4375 



0.4 

.6803 



0.16 

.4903 

.87461 


.11501 

-.15731 




-.08711 




-.08711 

.9353 

.34701 






-.4264 

.08561 

-.9840 

.37011 



HO 

-.3782 

.31871 

.3681 






.7055 

3.1925 



.3483 

2.8352 

.11181 






-.18981 

-.945411 



-.15701 

-.91221 

Auxiliary Matrix 












1.0169 

-0.3801 

0.2690 

-0.2258 



0 

0.6631 



0 

0.6631 

-.05241 

.05311 

-.06331 

.06611 




.05581 




.05581 

.9401 

.2221 

-1.7415 

2.0302 



1.7304 

3.0192 



0.6922 

1.9810 

-.05571 

-.04491 

.16581 

-.29791 



.34981 

.21821 



.13991 

.00031 

.6672 

.6713 

.9398 

-1.5922 



-.4285 

-1.0207 



.1278 

-.4643 

-.06841 

-.18521 

-.34871 

.08331 



-.06791 

.01511 



.08391 

.16681 

.2966 

.5880 

1,6566 

2.4870 



.2809 

1.2808 



.1458 

1.1456 

.09721 

-.02931 

-.34881 

-.75561 



.00901 

.00901 



-.01891 

-.01881 





k' 

1 

la 

= 0.5264 

1.5264 

k ' 

Ih 

» 0.3377 

1.3378 







.04851 

.04841 



.03201 

.03191 





k’ 

3 a 

« 1.1786 

2.1785 

k’ 

3 h 

* 1.0329 

2.0330 






.27791 

.27791 


.23491 

.23471 





k' 

5 a 

= 0.0195 

1.0193 

e' 

5 h 

* 0.3584 

1.3582 







-.07701 

-.07731 



.04171 

.04141 

NATIONAL 

ADVISORY 
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« 0.2809 

1.2808 
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K 

7 a 

k‘, 

7 h 

« 0.1458 

1.1456 







.00901 

.00901 



-.01891 

-.01881 
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TABLE XII.- VALUES OF w AT VARIOUS STATIONS ALONG THE SPAN 

(EXAMPLE I) 


y 



sin 

sin na( 

sin nd 



n 




toS 



n 

n 


n 


r 

1 

0.9165 

0.3095 + 0.02931 

0.4824 

+ 0.04451 

J 

' 

3 

-0.1099 

-0.1135 - 0.02581 

-0.1295 

- 0.03051 

0.4 < 


5 

-0.0936 

-0.0335 - 0.00391 

-0.0018 

+ 0.00721 



7 

0.1380 

0.0201 - 0.00261 

0.0388 

+ 0.00121 ' 



w; 

0.1826 - 0.00501 

0.3899 

+ 0.02241 



1 

0.800 

0.2702 + 0.02561 

0.4211 

+ 0.05881 

0.6< 


3 

0.1173 

0.1212 + 0.02761 

0.1382 

+ 0.03261 


5 

-0.1994 

-0.0715 - 0.00851 

-0.0039 

+ 0.01541 



7 

0.0295 

0.0043 - 0.00061 

0.0085 

+ 0.00031 



WJ 

0.3242 + 0.04451 

0.5637 

+ 0.08711 



1 

0.600 

0.2026 + 0.01921 

0.3158 

+ 0.02911 



3 

0.3120 

0.3223 + 0.07331 

0.3677 

+ 0.08671 

0.8 < 

1 

5 

-0.0151 

-0.0054 - 0.00061 

-0.0003 

+ 0.00121 

1 

L 

7 

-0.1398 

-0.0204 +0.00261 

-0.0595 

- 0.00151 



WJ 

0.4991 +0.09451 

0.6439 

+ 0.11571 


r 

1 

0.4359 

0.1472 +0.01391 

0.2295 

+ 0.02111 



3 

0.3255 

0.3562 +0.07651 

0.3836 

+ 0.09051 

0.9 S 


5 

0.1550 

0.0556 + 0.00651 

0.0030 

- 0.01191 



7 

-0.0022 

-0.0005 +0.00001 

-0.0006 

- 0.00001 



w; 

0.53S7 +0.09691 

0.6155 

+ 0.09971 
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TABLE XIII.- CALCULATION OF aj (j = a, h) FUNCTIONS 
(EXAMPLE I) 

PARABOLIC BENDING 


y 

LqS 


fj(y) 

"J 1 

iJ 

P 4- 

''J 

fj(y) 

Jo-iJ 

0 


0 

•H 

8 

1 

8 

0.5850 + 0.0309^ 

t 00 —00 1 

0.4 

0.1526 - 0.00301 

0.16 

0.1413 - 0.01881 

0.5850 + 0.03091 

0.0832 - 0.00661 

0.6 

0.5242 + 0.04451 

0.36 

-0.0994 + 0.12311 

0.5860 + 0.03091 

-0.0620 + 0.06891 

0.8 

0.4991 + 0.09451 

0.64 

-0.2202 + 0.14771 

0.5850 + 0.03091 

-0.1334 + 0.07961 

0.9 

0.5387 + 0.09691 

0.81 

-0.3349 + 0.11961 

0.5850 + 0.03091 

-0.1996 + 0.05961 

1.0 

0 

1.0 

- 1 

0.5850 + 0.03091 

-0.5850 - 0.0309i 

LINEAR TORSION 

0 


0 

CD — CXf i 

I 

0.5850 + 0.0309 

1 

1 00 — ool 

0.4 

0.3899 ^ 0.02241 

0.4 

-0.0252 + 0.05601 

0.5860 + 0.03091 

-0.0165 + 0.03201 

0.6 

0.5637 + 0.08711 

0.6 

-0.0605 + 0.14521 

0.5850 + 0.03091 

-0.0399 + 0.08311 

0.8 

0.6439 + 0.11571 

0.8 

-0.1951 + 0.14461 

0.5850 + 0.03091 

-0.1186 + 0.07861 

0.9 

0.6155 + 0.09971 

0.9 

-0.3161 + 0.11081 

0.5860 + 0.03091 

-0.1883 + 0.05511 

1.0 

0 

1.0 

- 1 

0.5850 + 0.03091 

-0.5850 - 0.03091 


TABLE XIV.- VALUES OF INTECSANDS FOR IBE IN' GRAPHICAL INTEGRATION 

(EXAMPLE I) 


z 


z* 


ajjZ* 

OhZ® 

0.4 

0.064 

0.0256 

0.0832 - 0.00661 

0.0021 - 0.00021 

0.0053 - 0.00041 

0.6 

0.216 

0.1296 

-0.0620 0.06891 

-0.0080 + 0.01031 

-0.0134 + 0.01491 

0.8 

0.512 

0.4096 

-0.1334 0.07961 

-0.0546 + 0.03261 

-0.0683 + 0.04081 

0.9 

0.729 

0.6561 

-0.1996 + 0.05961 

-0.1310 + 0.03911 

-0.1455 - 0.03091 

1.0 

1.0 

1.0 

-0.5850 - 0.03091 

-0.5850 - 0.03091 

-0.5850 + 0.03091 

z 

a 

z 

3 

z 

a 

3 

a z 

a 

a z 




a 

a 

a 

0.4 

0.16 

0.064 

-0.0165 + 0.0320i 

-0.0011 + 0.00201 

-0.0026 + 0.00511 

0.6 

0.36 

0.216 

-0.0399 + 0. 08311 

-0.0086 + 0.01791 

-0.0144 + 0.02991 

0.8 

0.64 

0.512 

-0.0016 + 0.07861 

-0.0607 + 0.04021 

-0.0759 + 0.05031 

0;9 

0.81 

0.729 

-0.1883 + 0.05511 

-0.1373 + 0.04021 

-0.1525 + 0.04461 

1.0 

1.0 

1.0 

-0.5850 - 0.03091 

-0.5850 - 0.03091 

-0.5850 - 0.03091 
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TABLE XV” 

AERODYNAMIC CORRECnONS AND CALCULATED DETERMINANT TERMS CORRECTED FOR FINITE-SPAN EFFECTS (EXAMPLE I) 

[ " • ( ^ ) ( 1 + 18 ) 6 = 8 „ - 8 ^] 
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TABLE r/I 


CHARACTERISTICS ALONG SPAl'I OF TAIL SURFACES USED IN EXAMPLE II 


y 

Horizontal 

tail 

1 

Tertlcal tall 1 


b ft ; 

1 

1 b ft 
1 

o 1 

! a 
1 

i 

i c 

j 

1 e 

i f 

0 

4 

1 

1 

4.64 

^ i 

1 

-1.909 

! -0.032 

j 

0.246 

I 

1 

i 1 
! 

.3 

— 

— 

3.95 I 

] 

1 .851 

1 

- 2.413 

1 

1 


1 .828 

.4 

3.44 

0.860 

j 

i 3.69 

j 

.797 

-2.634 

-.230 

.098 

1 .835 

i 

.513 



3.33 

.718 

1 

-2.970 

i 

1 


i .877 

.6 

3.08 

.770 

3. 00 

.648 

-3.300 i 

( 

-.540 

i 

.060 

1 

.3 

2.62 

.655 

2.18 

.471 i 

! -4.587 

-.404 

.037 

1 

1 

1 .® 

0 

0 

0 

0 1 
1 

— 

— 

0 1 

1 


TABLE XVII 


GENERAL CHARACTERISTICS OF TAIL SURFACES USED IN EXAMPLE II 



Horizontal tail 

i Vertical tail 

s 

2.525 ! 

1 

1 

i 0 1 

1 1.800 

i 


4 ft 

4.64 ft ! 


3.33 ft 

3.33 ft 

nb 

— 

-10.08 ft 

n 

i 

-■. 40 (lag) 

1 
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TABLE mil EVALUATION OF INTEJGRALS EQR EXAMPLE II 


^ ) U + Ig) where g -= = gj,] 


rsterminant 

term 

Structural peorb 
(ft) 

Two-dimens i onal 
aerodynamic part 

(ft) 

Detorminant term 
two-dimensional 
aerodynamic part 
(ft) 

' (k^ = 0.686 k^^ = 0.795) 

A 

67.27 [l - 1.988 n) 

2.1228 - lU. 66481 

69.393 - 14.66481 - 133.73 n 

B 

100.1 

-9.0606 - 24.933k! 

91.039 - 24.93341 

C 

1.731 

-11.2670 - 3.69931 

-9i5360 - 3.69931 

D 

100.1 

4.4524 - 17.83191 

104.552 - 17.83191 

E 

1891 [1 - n ] 

-27.7895 - 167.0931 

1863.2 - 167.0931 - 1891 n 

F 

23.862 

-67.1873 - 38.63701 

-93.3253 - 38.63761 

G 

1.737 

0.3202 - 0.17281 

2.0572 - 0.17281 

H 

23.96 

0.2137 - 4.7161 

24.1297 - 4.7161 

I 

5.99r 

-0.1613 - 3.25101 

5.8307 - 3.25101 


(k 

oH - 0.600 kov = 0.696) 


A 

67.27 [1 - 1.988 n] 

1.5917 - 17.39401 

68.862 - 17.39*^1 - 133.73 n 

B 

100.1 

-13.9834 - 28.86861 

86.117 - 28.86861 

C 

1.731 

-14.8607 - 3.57031 

-13.1297 - 3.57031 

D 

100.1 

3.7185 - 20.12901 

103,8185 - 20.12901 

E 

1891 ti - fl J 

-54.9543 - 184.6361 

1836 - 184.631 - 1891 n 

F 

23.862 

-88.6384 - 39.56451 

-64.7764 - 39.56451 

G 

1.737 

0.3108 - 0.20441 

2.0478 - 0.20441 

H 

23.916 

2.0841 - 5.31281 

26.0000 - 5.31281 

I 

5.992 

-0.4069 - 3.50801 

5.5851 - 3.50801 


(k 

OH = = 0 - 506 ) 


A 

67.27 [1 - 1.988 n] 

-0.6677 - 25.85681 

66.6023 - 25.8568 - 133.73 n 

B 

100.1 

-34.9401 - 39.19341 

65.1599 - 39.19341 

C 

1.731 

-30.2362 - 2.94831 

-28.5052 - 2.94831 

D 

100.1 

0.7482 - 29.57011 

100.8482 - 29.57011 

E 

1891 [ 1 - n ) 

-166.9957 - 255.16971 

1724 - 255.16971 - 1891 n 

F 

23.862 

-177.2376 - 43.31121 

-153.3756 - 43.3U21 

G 

1.737 

0.2795 - 0. '28971 

2.0165 - 0.28971 

H 

23.916 

0.7526 - 7.24m 

24.6686 - 7*. 24m 

I 

5.992 

• -i.io38 - 4.72781 

4.8882 - 4.72781 
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TABLE XIX.- SOLUTION OF SET OF EQUATIONS 

THE GROUT METHOD (EXAMPLE II ) 


(a) 


Given Matrix 



Oh 

Check 


O.I4.7O9 

-0.11941 

- 0 . 4 i ;59 

0.15661 

0.2507 

-0.09521 

0.0591 

0.55181 

0.5348 

0.49581 



•-0.2670 

0.09781 

- 0.2849 

0.11111 

0.0768 

0.77091 

0.1808 

0.79541 


0.6347 

-0.17601 

0.2479 

-0.08621 

-0.2657 

0.11271 

0.0996 

0.92451 

0.7165 

0.77501 

Supplementary Matrix 




1.9955 

o. 5059 i 


-0.4459 

0.15661 

0.2507 

-0.09521 

0.0591 

0.55181 

0.5548 

0.49581 

2.^5991 

o.9625i 



-0.6350 

0.25081 

-0.0260 
-0. 00061 

-0.3065 

0.11241 

0.9675 

0.5l].i6i 




0.0649 

0.18091 

0.9845 

-0.14591 

Auxiliary Matrix 






0.4709 
-0. 11941 

-0.9689 

0.08691 

0.5474 

-0.05911 

-0.1612 

1.15091 

0.4172 

1.15861 


0.6599 

-o.i 844 i 

0.5525 

-0.15791 

-1.7985 

0.00781 

-0.0634 

-0.02651 

-0.86211 

-0.01851 


0.6547 

-0.1760 

0.8476 

-0.51591 

0.9192 

-0.52461 

0.0010 

0.19721 

1.0015 

0.19701 




^2h 

-0.2051 1 
1 . 54621 

0.7970 

1.54601 
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^6h 

-0.0601 

0.52821 

0.0010 

0.19721 

0.9405 

0.32801 

1.0015 

0.19701 
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TABLE XX,- VALUES OF Q FOR EXAMPLE II 


y 

1 >o® 

n 

sin n<f 
n 

sin n^ 

1 


2 

0.5667 

-0.0745+0.49571 

0.4^ 


k 

-0.2i;93 

0.0150-0.08181 


L 

6 

+0.1152 

+0.0001+0.02251 




-0.0594+0.45421 

I 

r 

2 

o.ij.800 

-•0.0975+0.64621 

o.6<j 


k 

-0.15104. 

0.0081-0.04411 


L 

6 

-0.1098 

-0.0001-0.02171 




-0.0895+0.58041 

I 

r 

2 

0.4800 

-0.0975+0.64621 

0.8^ 

1 

k 

0.1544 

-0.0081+0.04411 

1 

L 

6 

-0.1098 

-0.0001-0.02171 



Oh 

-0.1057+0.66861 


TABLE XXI,- AERODYNAMIC CORRECTIONS AND CALCULATED DETERMINANT TERNB 
CORRECTED FOR FINITE-SPAN EFFECTS (EXAMPLE II) 



koH = 0.436 

kbV 

= 0,506 

Term 

Det. term with two- 
dim. aero. part. 

Finite “Span 
corrections 

Det. term corrected for 
finite-span effects 

A 

66.6025-25.85681 

-155.75 n 

5 .0885+6 . 98841 

71.6908-18.86841 
-155.75 n 

B 

65.1599-59.19541 

11.5520+2.86801 

76.7119-56.52541 

C 

-28.5052-2.94851 

6.5595-5.25601 

-21.9459-6.18451 

D 

100.8482-29.57011 

6.5861+9.04921 

107.2545-20.52091 

E 

1724-255,16971 
-1891 n 

65.4194+16.16561 

1789.4-259.0041 
-1891 a 

F 

-155.5756-45.51121 

54.6550-17.11151 

-118.7415-60.42251 

G 

2.0165-0.28971 

0.0572+0.09461 

2.0757-0.19511 

H 

24.6686-7.24111 

-4.7027+1.75^^-61 

19.9659-5.48651 

I 

4.8882-4.72781 

0.5885-0.20121 

5.2767-4.92901 
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TABLE XXI.— Concluded 


^oH " 0*600 k^y = 0.696 


A 

68.8617-17.59401 
-155.75 n 

2.8478+5.98221 

71. 7095-15. 41181 
-155.75 n 

B 

86.117-28.86861 

5.8026+2.49851 

91 .9196-26.57051 

C 

-15.1297-3.57051 

2.8707-1.10771 

-10.2590-4.67801 

D 

105.819-20.12901 

5.0659+5.64251 

106.8829-14.48671 

E 

1856-18I4-.6511 
-1891 0 

29.5925+15.01881 

1865. 6-169.6171 
-1891 n 

F 

- 614 -. 7764-59 .5^51 

14.5951-5.28151 

-50.1855-44.84601 

G 

2.014.78-0.20441 

0.0264+0.05611 

2.0742-0.14851 

H 

; 26.0001-5.51281 

-1.9255+0.52501 

24.0768-4.98781 

I 

5.5851-5.50801 

0.1624-0.04181 

5.7475-5.54981 



0.686 k^^ = 

0.795 1 

A 

69.5926-14.66481 

-155.75 n 

1.9815+5.06751 

71.5741-11.59751 
-155.75 fj 

B 

91.05-94-24.95541 

4.0857+2.19551 

95.1251-22.75791 

C 

-9.5560-5.69951 

1.9556-0.59261 

-7.6024-4.29191 

D 

104.5524-17.^191 

2.2890+4.26181 

106.8414-15.57011 

E 

1865.2-167.0951 

-1891 n 

20.5217+15.12001 

1885. 7-155 . 9751 

-1891 n 

P 

-45.5255-58.65701 

9.8602-2.50801 

55.4651-41.14501 

G 

2.0572-0.17281 

0.0182+0.04161 

2.0754-0.15121 

H 

26.055-4.7161 

-1.2895+0.09551 

24.7657-4.62251 

I 

5.8507-5.25101 

0.1060-0.01721 

5.9567-5.26821 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS. 




Figure 1.- (a), Plan form of wing, showing location of elastic axis and midchord line. 

(h), Wing section with aileron and tab, showing main parameters in units of 
local semichord. 
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Fig. 2a 
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(a) Cos 0=1; n = 1, 2, 3, 4, 5, 6, 7, 8, 9. 


Figure 2 (a to h.) 
imaginary parts; 


- Function Sn(0, ICq®) against koS for various 
values of n and of cos 0. (Dash lines are lor 
)lid lines, for real parts.) 
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Fig. 2b 



(b) Coe 0 = 0.8} n = 1, 3, 5, 7, 9. 
Figure 2 .- Continued. 


S.(»>«s) 


Fig. 2c 
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Figure 2.- Continued. 
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Fig. 2d 



Figure 8.- Continued. 




Fig. 2e 
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Figure 2.- Continued. 
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Fig. 2f 



koS 

(f) Cos 0 = 0.4j n = 1, 3, 5, 7, 9. 


Figure 2.- Continued. 


Fig. 2g 
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kos 

(g) Cos 0 = 0.?; n = 1, 3, 7, 9. 

Figure 2 .- Continued. 
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Fig. 2h 



(a) Cos = Oj n = 1, 3, 5, 7, 9, 


Figure 2.- Concluded. 



Fig. 3 
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/« = 3 


/R = 6 
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figure 3.- Lift distributions for uniform translation of rectangular 
v^ings of aspect ratios 3 and 6 for ko = 0.167, 0.333, 0.667. 
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Figure 4.- Lift distributions for uniform pitching of rec- 
tangular wings of aspect ratios 3 2Uid 6 for Icq = 
0, 0.167, 0.335, 0.667. 
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5.- Lift distributions for linear symmetrical torsion 
of rectangular wings of aspect ratios 3 and 6 for 
0.167, 0.333, 0.667. 
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figure 6.- 
ko = 0.167 


Lift distributions for linear symmetrical bending 
of rectangular wings of aspect ratios 3 and 6 for 
0.333, 0.667. 
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7.- Lift distributions for parabolic symmetrical bend- 
ing of rectangular wings of aspect ratios -3 and 6 
= 0.167, 0.333, 0.667. 
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Figure 8.- Lift distributions for linear axiti symmetrical torsion of a 

rectangular wing of aspect ratio 3 with = 0, 0.16V, 0.333, 
0.667 (left column), and for rigid rolling (right column) of a rectangular 
wing of aspect ratio 3 with = 0.167, 0.333, 0.667, 1.333.' 
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Figure 10.- 
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Lift distributions 
of elliptical wings 
0.424, 0.847. 


for linear symmetrical bending 
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figure 11.- Lift distributions for parabolic symmetrical bend- 
ing of elliptical wings of aspect ratios 3 and 6 
with ko = 0.212, 0.424, 0.847. 


NACA TN No. 1195 


Fig. 12 


/R=3 ^ = ^ 



Figure 12.- Lift distributions for linear antisynunetrical torsion 
of an elliptical wing of aspect ratio 3 with k© * 0, 
0.212, 0.424, 0.847 (left column; and for rigid rolling of an 
elliptical wing of aspect ratio 3 with ko = 0.212, 0.424, 0.847, 
1.695 (right colximn) . 
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Figure 13.- Lift distributions and hinge-raoment distributions 
for full-span aileron v/ith uniform spanwise de- 
flection for an elliptical wing of aspect ratio 3 (e = C, 1 = 
0.2b/bo) ko = 0, 0.424, 0.847. 
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Figure 14,- Lift and hinge-moment distributions for an ellipti- 
cal wing of aspect ratio 6 with half-span ailerons 
of approximately one-quarter the wing chord and with ko = 0, 
0.424, 0.847 (X = O.lb/bo). 
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^''igure 16.- Rolling-moment coefficients for an elliptical wing of aspect 
ratio 3 in rigid rolling and for an elliptical wing of aspect 
latio 3 in linear antisymmetrical torsion plotted as vector diagrams and 
as functions of ko . 
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Figure 17.- Lift coefficients for an elliptical wing of aspect ratio 3 
in uniform pitching and for an elliptical wing in uniform 
vertical translation plotted as vector diagrams and as functions of ko- 
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Ficure 18.- Lift distributions for uniform translation of an elliptical 
wing of aspect ratio 3 for k® = 0.212, 0.424, O.S**?, 1.695. 
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Figure 19.- Lift distributions for uniform pitching of an elliptical 

wing of aspect ratio 3 for ko = 0, 0.212, 0.424, 0’.847, 1.695. 
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Figure 22 (a to d).- Unifor/n pitching and translation of rigid elliptical wing (one-point 
collocation). Comparison of functions C , C + cr with functions F and 

Q of reference 5. 
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Figure 20.- First-torsion, first-bending, and second-bending modes for a 
rectangular wing of aspect ratio 13.5 and lift distributions 
for these modes for k© = 0.0965. 



Figure 21.- Lift distributions for linear eintisymmetrical torsion of an 
elliptical wing of aspect ratio 3 for Icq = 1.695 and 2.542. 
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Fig. 22c ,d 




Figure 22.- Concluded. 
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Figure 23.- Damping coefficient against true airspeed for a 
rectangular wing of aspect ratio 6. b© = 10 
inches; g = 0.070, (Analysis by method of reference 4.) 
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Figure 24.- Concluded, (d) For AE: C 
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Figure 25.- Horizontal tall used in example II. 
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Figure 26.- Fin-rudder combination used in example II. 
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Figure 27.- Damping coefficient g against true airspeed for 

the tail-flutter problem of example II; g = 0.038 
Note« Flutter determinants solved by method of reference 4. 


